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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 
RESEARCH  MEMORANDUM 

DAMPING  IN  PITCH  OF  LOW-ASPECT-RATIO  WINGS 
AT  SUBSONIC  AND  SUPERSONIC  SPEEDS 
By  Murray  Tobak 

SUMMARY 


The  concept  of  indie ial  functions  is  applied  to  the  analysis  of  the 
aerodynamic  phenomena  associated  with  the  short-period  pitching  mode  of 
wings  in  subsonic  and  supersonic  flight.  Simple  physical  relationships 
are  pointed  out  and  are  used  to  study  the  effect  on  the  rotary-damp ing- 
moment  coefficient  of  changes  in  center-of -gravity  position,  Mach  number, 
aspect  ratio,  plan  form,  frequency,  and  thickness.  Qualitative  conclu- 
sions are  drawn  from  the  results  of  this  investigation  and  are  compared 
with  the  results  of  experiments  for  a series  of  low-aspect-ratio  wing- 
body  combinations  having  triangular,  swept,  and  unswept  wing  plan  forms . 

Results  of  the  experimental  investigation,  which  were  obtained  by  a 
single-degree-of-freedom  free-oscillation  technique  over  the  Mach  number 
ranges  0.6  to  0.9  and  1.2  to  1.9>  were  in  good  agreement  with  the  results 
of  theoretical  computations.  The  predictions  of  ranges  of  supersonic 
Mach  numl er  and  center-of-gravity  positions  over  which  dynamic  instabil- 
ity may  be  expected,  of  the  beneficial  effect  on  the  damping  moment  of  a 
reduction  in  aspect  ratio,  and  of  only  a small  effect  of  thickness  on  the 
damping  moment  were  confirmed  by  the  experimental  results. 

The  occurrence  at  high  subsonic  Mach  numbers  of  small-amplitude 
self-sustained  pitching  oscillations  is  noted,  and  a hypothesis  is 
advanced  for  its  explanation. 


INTRODUCTION 


In  the  classical  study  of  the  longitudinal  motion  of  an  aircraft, 
it  is  usually  found  that  the  motion  resulting  from  a small  equilibrium- 
destroying  disturbance  consists  of  two  modes:  one,  a lightly  damped,  low- 
frequency  motion  at  essentially  constant  attitude,  called  the  phugoid 
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oscillation;  the  other,  a rotary-pitching  and  plunging  oscillation  of 
high  frequency  (relative  to  the  phugoid  frequency)  called  the  short- 
period  oscillation.  The  phugoid  oscillation  has  generally  been  described 
as  resulting  from  a slow  interchange  of  potential  and  kinetic  energy  as 
the  aircraft  experiences  periodic  variations  in  airspeed  and  altitude. 

The  character  of  the  phugoid  motion  as  influenced  by  airspeed,  altitude, 
and  aircraft  geometry  has  been  well  understood  for  some  time  (see,  e.g., 
ref.  l) . The  short-period  motion,  on  the  other  hand,  having  in  the  past 

been  found  to  be  highly  damped  and  of  short  duration,  has  been  the  cause 

of  no  concern.  Its  characteristics  therefore  have  not  been  as  fully 
investigated  as  those  of  the  phugoid  oscillation.  With  the  advent  of 
flight  at  speeds  approaching  and  exceeding  the  speed  of  sound,  however, 
the  loss  of  rotary  damping  occurring  in  practically  all  aircraft  at 
speeds  near  the  sonic  speed  has  caused  renewed  interest  in  the  short- 
period  pitching  mode.  Unlike  the  easily  controlled  phugoid  oscillation, 
the  deterioration  of  damping  in  the  short-period  mode  is  of  serious  con- 
cern to  the  pilot,  since  the  period  of  the  oscillation  can  be  of  the 
same  order  of  magnitude  as  the  pilot's  reaction  time.  The  oscillation 

may  therefore  be  difficult  or  even  impossible  for  the  pilot  to  control 

manually.  Furthermore,  the  additional  load  imposed  upon  the  airframe 
due  to  a rapid  growth  of  the  amplitude  of  a negatively  damped  oscillation 
makes  possible  the  occurrence  of  structural  failure.  It  is  therefore  of 
considerable  interest  to  obtain  an  understanding  of  the  nature  of  the 
short-period  mode,  parallel  to  that  which  has  been  gained  of  the  phugoid 
mode . 


One  means  of  viewing  the  aerodynamic  phenomena  occurring  during  the 
short-period  oscillation  from  a fundamental  standpoint  is  through  appli- 
cation of  the  concept  of  indicial  functions.  In  this  approach,  the  var- 
iations with  time  of  the  aircraft  angle  of  attack  and  angular  velocity 
during  the  oscillation  are  replaced  by  a large  number  of  small  instanta- 
neous or  step  changes.  The  transient  aerodynamic  reactions  to  these  step 
changes  are  termed  "indicial  functions,"  and  have  been  calculated  theo- 
retically for  several  classes  of  wings  (refs.  2 to  6) . By  suitable 
superposition  of  these  results  (refs.  7 to  9)>  the  aerodynamic  forces  and 
moments  caused  by  the  given  maneuver  can  be  studied.  It  will  be  the 
purpose  of  this  report  to  make  such  a study  for  the  simplified  case  of  an 
aircraft  performing  single-degree-of-freedom  rotary  oscillations.  For 
this  maneuver,  which  corresponds  to  the  short-period  oscillation  when  the 
plunging  velocity  of  the  aircraft  is  zero,  the  use  of  simple  physical 
relationships  associated  with  the  indicial  function  concept  enables 
qualitative  studies  to  be  made  of  the  separate  effects  on  the  aerodynamic 
forces  and  moments  of  changes  in  Mach  number,  aspect  ratio,  plan  form, 
frequency,  and  thickness.  Results  of  this  investigation  are  then  com- 
pared with  the  results  of  experiments  with  a group  of  low-aspect-ratio 
wing-body  combinations.  The  tests  were  conducted  in  the  Ames  6-  by  6-foot 
supersonic  wind  tunnel  and  were  similar  in  technique  to  those  reported  in 
reference  10. 
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NOTATION 


g 

A aspect  ratio,  b /S 

Cl  lift  coefficient,  lift/qgS 

Cm  pitching -moment  coefficient,  pitching  moment/c^Sc 

P-P0 

Cp  pressure  coefficient,  — - 

I moment  of  inertia,  slug-ft2 

Mq  free-stream  Mach  number,  V0/a0 

R Reynolds  number,  based  on  wing  mean  aerodynamic  chord 

S wing  area,  including  portion  enclosed  by  body,  sq  ft 

V0  flight  speed,  ft/ sec 

a0  speed  of  sound  in  free  stream,  ft/sec 

b wing  span,  ft 

c0  wing  root  chord,  ft 

2 b/2  2 

c wing  mean  aerodynamic  chord,  — / (local  chord)  dy 

D O 

e base  of  natural  logarithms 

i /r 

k reduced  frequency  parameter,  wc/2V0 

Ap  local  loading  at  plan-form  surface, 

pressure  lower  surface  - pressure  upper  surface 

% 

q angular  velocity  due  to  pitching,  radians/sec 

q0  free-stream  dynamic  pressure,  ^ pQV02,  lb/sq  ft 
t time,  sec 

ta  time  required  following  an  instantaneous  change  in  angle  of  attack 
or  angular  velocity  for  the  transient  lift  or  moment  to  attain 
steady  state,  sec 
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x,y,z  rectangular  coordinates 

xa.c.  distance  from  leading  edge  of  M.A.C.  to  aerodynamic  center,  ft 
x0  distance  from  leading  edge  of  M.A.C.  to  axis  of  rotation,  ft 
AXq  xa.c.  - xc 

a angle  of  attack  of  wing  center  line  with  respect  to  free-stream 

direction  (sketch  a) 

3 yiMo2-ii 

7 ratio  of  specific  heat  at  constant  pressure  to  that  at  constant 
volume 

5 airfoil-thickness  ratio,  maximum  thickness/chord 

0 angle  of  wing  center  line  with  respect  to  horizontal  axis 
(sketch  a) 

V acute  angle  between  wing  plane  of  symmetry  and  trailing  edge 
(sketch  p) 

p0  free-stream  density,  slugs/cu  ft 

to  angular  frequency  of  oscillation,  radians/sec 

oa  distance  traveled,  measured  in  half  M.A.C.  lengths,  in  the  time 
interval  ta,  2V0ta/5 

When  a,  a,  and  q are  used  as  subscripts,  a nondimens ional  deriv- 
ative is  indicated,  and  this  derivative  is  evaluated  as  the  independent 
variable  (a, a,  or  q)  approaches  zero.  For  example. 
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THEORY 

Application  of  Indicial  Functions  to  the  Aerodynamic 
Theory  of  Unsteady  Flows 


One  of  the  most  useful  tools  in  the  study  of  unsteady  flows  is  the 
concept  of  indicial  aerodynamic  functions  which  may  be  defined  briefly 
as  the  aerodynamic  response  of  the  airfoil  as  a function  of  time  to  an 
instantaneous  change  in  one  of  the  conditions  determining  the  aerody- 
namic properties  of  the  airfoil  in  a steady  flow.  Theoretical  aerody- 
namic indicial  functions  were  first  derived  by  Wagner  (ref.  2)  for  the 
two-dimensional  wing  in  incompressible  flow.  More  recently,  these 
results  have  been  extended  by  Heaslet  and  Lomax  to  cover  the  compressible 
case  for  both  subsonic  and  supersonic  speeds  (ref.  4).  In  addition, 
theoretical  indicial  functions  have  now  been  obtained  for  both  wide  and 
slender  triangular  wings  and  rectangular  wings,  all  for  supersonic  speeds 
(refs.  4 to  6). 

The  indicial  function  derives  its  usefulness  primarily  through  the 
ease  with  which  it  lends  itself  to  the  powerful  and  well-established 
methods  of  the  operational  calculus  (refs. 7 to  9)*  With  the  use  of  these 
methods,  the  aerodynamic  forces  and  moments  caused  by  arbitrary  motions 
of  the  airframe  can  be  studied  from  a fundamental  standpoint.  Because  of 
the  wide  range  of  applicability  of  this  means  of  approach  in  unsteady  flow 
analyses,  a considerable  portion  of  the  succeeding  discussibn  is  devoted 
to  the  fundamentals  involved. 

Definition  of  coordinate  system.-  In  the  succeeding  analysis  the 
stability  system  of  axes  is  used.  The  origin  of  the  coordinate  system 
is  placed  in  the  airfoil  so  that  the  y axis  which  is  perpendicular  to 
the  vertical  plane  of  symmetry  is  coincident  with  the  axis  of  rotation  of 
the  airfoil;  the  positive  branch  of  the  x axis  is  pointed  in  the  direc- 
tion of  flight;  and  the  z axis  lies  in  the  vertical  plane  of  symmetry, 
positive  downward.  The  angle  of  attack  a is  measured  as  the  angle 
between  the  chord  plane  of  the  airfoil  and  the  xy  plane,  and  is  shown 
as  positive  in  sketch  (a).  The  angle  of  pitch  6 is  the  angle  between 
the  chord  plane  of  the  airfoil 
and  the  horizontal  plane  (an 
arbitrary  reference)  and  is  also 
shown  positive  in  sketch  (a). 

Forces  are  measured  as  positive 
upward,  whereas  pitching  moments 
are  positive  when  tending  to 
increase  the  angle  of  pitch  in 
the  positive  direction.  When 
the  airspeed  V0  is  constant 
which  corresponds  to  the 
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condition  under  study,  the  translatory  and  angular  motions  of  the  air- 
foil with  respect  to  any  system  of  coordinates  are  defined  if  the  time 
histories  of  the  angle  of  attack  a and  the  angle  of  pitch  0 and 
their  derivatives  are  known.  For  purposes  of  clarity,  two  different 
harmonic  motions  of  the  aircraft  are  shown,  illustrating  the  difference 
between  a flight  path  which  involves  a constant  angle  of  attack  and  a 
varying  angle  of  pitch  (sketch  (b))  and  one  which  involves  a constant 
angle  of  pitch  and  a varying  angle  of  attack  (sketch  (c)). 


of  attack  • O angle  of  attack  • a 

Sketch  (b)  Sketch  (c) 


Now  consider  the  case  of  a wing  executing  harmonic  rotary  oscilla- 
tions about  the  y axis  while  the  origin  of  the  coordinate  system 
traverses  a level  path  at  constant  velocity  VQ.  This  case  corresponds 
to  that  of  a wind-tunnel  model  mounted  to  permit  single-degree-of- 
freedom  rotary  oscillations,  or  to  the  short-period  mode  of  an  aircraft 
in  flight  when  the  plunging  velocity  of  the  center  of  gravity  is  zero. 
Here  a and  9 are  equal,  so  that  the  maneuver  is  defined  by  one  vari- 
able, the  time  history  of  either  a or  0.  Let  the  angle  of  attack 
be  a and  the  angular  velocity  be  q (q  s d0/dt  = da/dt).  At  any 
instant,  the  normal  velocity  at  any  point  on  the  airfoil  surface  is  com- 
posed of  two  parts,  one  due  to  the  instantaneous  angle  of  attack  aVQ, 
the  other  due  to  the  angular  velocity  at  the  same  instant  -qx  (see 
sketch  (d)).  These  are  two  of  the  instantaneous  boundary  conditions  of 
the  unsteady  flow. 


normal  velocity 
kve  to  angle  of  attack 
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Solutions  for  the  aerodynamic  forces  and  moments  which  correspond 
to  these  boundary  conditions  may  be  derived  by  a number  of  methods 
involving  various  degrees  of  approximation.  In  succeeding  sections, 
the  use  of  the  concept  of  indicial  functions  and  the  principle  of  super- 
position for  this  purpose  will  be  illustrated  and  compared  with  other 
current  widely  used  methods. 

Concept  of  indicial  functions.-  In  order  to  illustrate  this  concept, 
assume  that  the  airfoil  under  consideration  has  been  flying  a level  path 
at  zero  angle  of  attack.  At  some  time,  which  is  designated  time  zero, 
the  wing  is  caused  to  attain  simultaneously  a constant  angle  of  attack  a 
and  angular  velocity  q.  The  normal  velocity  of  the  flow  next  to  the 
surface  of  the  airfoil  therefore  changes  discontinuously  from  zero  to  a 
pattern  that  is  constant  with  time  and  identical  in  shape  to  the  pattern 
shown  previously  in  sketch  (d) . The  lift  and  pitching  moment  that 
result  are  of  a transient  character  and  attain  their  steady-state  values 
corresponding  to  these  new  boundary  conditions  only  after  a significant 
interval  of  time  has  passed.  It  should  be  noted  there  exists  an  essen- 
tial difference  between  the  length  of  this  time  interval  at  subsonic  and 
supersonic  speeds.  At  supersonic  speeds,  the  vorticity  shed  into  the 
airfoil  wake  cannot  influence  the  flow  about  the  airfoil  but  at  subsonic 
speeds  this  influence  exists  for  all  time.  The  result  is  that  the  lift 
and  moment  reach  steady-state  values  in  a finite  time  at  supersonic 
speeds  but  approach  these  values  asymptotically  at  subsonic  speeds.  In 
either  case,  however,  the  time  responses  in  lift  and  moment  to  the  step 
changes  in  a and  q are  termed  indicial  functions.  Sketch  (e)  illus- 
trates typical  subsonic  and  supersonic  indicial  lift  responses  to  a step 
change  in  the  angle  of  attack. 


Sketch  (e) 
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It  is  obvious  that  the  time  history  of  the  wing  motion  during  a 
short-period  oscillation  may  be  broken  down  into  an  infinite  number  of 
infinitesimally  small  step  changes  in  the  angle  of  attack  and  step 
changes  in  the  angular  velocity.  The  summation  of  the  indicial  lift  and 
moment  for  these  steps  then  yields  the  total  lift  and  moment  at  any  pre- 
scribed time.  In  sketch  (f),  the  mechanics  of  the  procedure  are  illus- 
trated for  an  arbitrary  angle-of -attack  variation.  Here,  the  given 

angle-of-attack  variation  is 
replaced  by  a number  of  small 
step  changes.  Within  each  step 
the  corresponding  response  in 
lift  is  shown  plotted  for  con- 
venience. It  is  then  apparent 
that  the  total  lift  at  time  t 
is  equal  to  the  sum  of  the 
increments  of  lift  in  each  step 
at  time  t.  As  indicated  by 
the  leaders,  however,  it  is 
clear  that  the  increments  of 
lift  for  the  various  steps  at 
time  t are  equivalent  to 
increments  in  the  first  step 


at  time  t 


Alterna- 


tively, then,  the  total  lift 
at  time  t can  be  written  as1 


Sketch  (f) 

t 

Cx,(t)  = CL^tMo)  + ^ CLaCt-ti)  (t^Atj^  (l) 

o 

After  a transformation  of  variables,  t - tj_  = T , and  letting  the  incre- 
ment of  time  approach  zero,  equation  (l)  can  be  rewritten  in  a form  of 
Duhamel's  integral  (see,  e.g.,  ref.  9) 


CL(t)  = “ J*  CLq,( *r) a( t— t) dT  (2) 

^■Here,  and  in  the  remainder  of  this  report,  the  use  of  parentheses  is 
reserved  solely  for  the  indication  of  functional  dependence.  Thus, 
for  example,  in  equation  (l)  the  term  Cj^Ct-ti)  is  interpreted  as 

the  value  of  Ct_  at  time  t - ti,  whereas  — (tj.)  is  the  value 
XJ«-  At 

of  Ax/At  at  time  All  other  enclosures  indicate  algebraic 

expressions  in  the  usual  sense. 
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A similar  procedure  is  carried  out  for  the  angular  velocity 
variation,  whereupon  the  total  lift  coefficient  at  the  prescribed 
time  t becomes 

cL(t)  - St  £ °tu(T)a(t-T)dT  + St  £ %(T)  w0  i(t-T)dT  <3) 

It  should  be  pointed  out  that  in  this  form  equation  ( 3)  is  applicable 
to  the  analysis  of  arbitrary  motions,  the  only  restriction  being  that 
the  flight  speed  is  constant.  In  the  following  sections,  however,  the 
application  of  equation  (3)  is  restricted  to  harmonic  motions  having  a 
single  degree  of  freedom.  The  reasons  for  this  restriction  are  two- 
fold: first,  the  motions  of  a statically  stable  aircraft  in  response  to 
a disturbance  are  most  generally  of  a harmonic  nature;  and  second,  such 
a restriction  permits  an  assessment  of  the  influence  of  the  time  rate 
of  the  airfoil  motions  on  the  aerodynamic  forces  and  moments. 

Application  of  indicial  functions  to  harmonic  pitching  oscillations 
Consider  first  a pure  sinusoidal  pitching  oscillation,  the  angle  of 
attack  being  zero  throughout  the  motion.  The  flight  path  for  such  a 
motion  has  been  illustrated  in  sketch  (b).  In  this  case,  the  angle  of 
pitch  is  given  by 

e(t)  = e0eiwt 

where  Qq  is  the  maximum  amplitude  of  oscillation  and  w is  the 
angular  frequency.  The  angular  velocity  is,  of  course,  q(t)  = b = 
iu)0oeiwt  = iw0(t).  Inserting  the  value  for  q(t)  in  equation  (3)  and 
performing  the  indicated  operations,  there  results 


°L(t)  * - | ^ «<*>  £ CLq(T)€“dwTdT  + «<>%(*)  (») 
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Note  In  sketch  (g)  that  Cl^(t)  is  equal  to  CL^(t)-F2(T),  and  that  for 
subsonic  speeds  F2(t)  approaches  zero  as  t approaches  t. 


Sketch  (g) 

Replacing  Cl^t)  in  equation  (4)  by  this  equality, 


Cl( t ) Iwc 


2- 

0)C 


e(t) 


= Cr  (t)  + / F2(t) 

2v  "a  ow  JU  ^ 


-iWT 


2Vr 


dT 


(5) 


For  subsonic  speeds,  let  t approach  infinity.  With  this  sub- 
stitution, equation  (5)  thereby  represents  the  lift  coefficient  due  to 
the  harmonic  pitching  motion  after  the  transient  loading  subsequent  to 
the  start  of  the  motion  has  reached  steady  state.  Then  separating 
equation  (5)  into  components  in-phase  (real  part)  and  out-of-phase 
(imaginary  part)  with  0,  there  is  obtained 


Cl  w2c  - " _ , N , 

— = — — / Fj,(t)  COS  WTdT 

0 2V0  Jo  2' 


CL 

iwc 0/2Vo 


w J F2(t)  sin  wTdT 


(6) 


Introduce  the  nondimens ional  parameters. 


9 


T 


number  of  half  M.A.C.  lengths  traveled  in  time  t 


k 


2V0 


reduced  frequency 
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In  terms  of  these  parameters,  equation  (6)  becomes,  for  Mq<1, 

— = k2  (m  F2(q>)  cos  k<pd<p 
9 Jo 


cl 


— ; = CLq(»)  - k / F2(<P)  sin  k<pd<p 

qc/2V0  ^ Jo  J 


(7) 


At  supersonic  speeds,  equations  (7)  may  be  simplified  somewhat 
since  the  build-up  in  lift  is  completed  in  a finite  number  of  half 
M.A.C.  lengths  of  travel  oa.  In  equations  (7),  therefore,  the  upper 
limits  of  the  integrals  may  be  replaced  by  cra,  since  beyond  that  point 
F2(<p)  is  zero. 


Ct  ^ Q, 

— = k 2 / F2(<p)  cos  k<pd<p 

9 J0 


Cl 

qS/2V0 


* cLq(°a) 


<p  > a 

M0>1 


sin  k<pd<p 


(8) 


Thus,  it  appears  from  equations  (7)  and  (8)  that  there  are  both  in-phase 
and  out-of-phase  lift  forces  associated  with  the  harmonic  pitching  oscil- 
lation. Notice,  however,  in  equations  (8)  that  if  the  trigonometric 
terms  are  expanded  and  the  reduced  frequency  is  required  to  be  small 
compared  to  unity  (corresponding  to  the  frequencies  encountered  in 
dynamic  stability  work)  terms  containing  second  and  higher  powers  of  k 
will  be  very  small  compared  to  first-order  terms.  Thus,  for  slow  fre- 
quencies, the  only  force  of  consequence  during  the  pitching  oscillation 
is  the  first  order  in  frequency  out-of -phase  lift  force,  qc/2V0  CL^(aa).2 


aThis  quantity  is,  of  course,  the  same  lift  force  in  phase  with  the 
pitching  velocity  which  would  occur  alone  had  the  wing  been  executing 
a steady  turn  (q  constant).  One  of  the  chief  advantages  of  the  indi- 
cial  response  method,  at  least  for  supersonic  speeds,  is  the  ease  with 
which  the  relative  importance  of  the  various  terms  contributing  to  the 
total  lift  and  moment  can  be  assessed  and  the  sources  of  the  Important 
contributions  identified. 
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2V0  9 ^ 


Sketch  (h) 


The  phase  relationships  for  the  har- 
monic pitching  oscillation  are  indi- 
cated in  sketch  (h) . It  is  evident 
from  sketch  (h)  that  the  total  lift 
leads  the  angle  of  pitch  by  nearly  90°. 

Application  of  indicial  functions 
to  harmonic  plunging  oscillations . - 
Next,  consider  a purely  sinusoidal 
variation  of  the  angle  of  attack,  the 
angle  of  pitch  being  zero  throughout 
the  motion.  The  flight  path  for  this 
motion  has  been  illustrated  in 
sketch  (c).  Here,  a equals  otQe^10^, 
where,  as  previously,  ccq  and  u>  are 
the  maximum  amplitude  and  angular  fre- 
quency, respectively.  Applying  equa- 
tion (3)  again. 


°L(t>  ma»&£  ^(Tje^-^dT 


= iwa(t)  /*  CLaMe'^dT  + (9) 

Now,  as  in  the  previous  example,  let  Cj^Ct)  = CjJa(t)  - F^t)  so  that 
equation  (9)  becomes 

CL(t)  = aftJCLjt)  - iU)a(t)  / * Fi(T)e"iWT  dT  (10) 


Again,  introduce  the  nondimens ional  parameters  and  k,  separate 
equation  (10)  into  its  real  and  imaginary  parts,  and  let  <p  approach 
infinity  for  subsonic  speeds  and  oa  for  supersonic  speeds.  There 
results 


~ = cLa(^)  - * £ Fx(<p)  sin  k<pd<p 


CL 

a5/2V0 


f Fi(^)  008  k<PecP 


X = »,  Mq<1  I 
\ * aa,  M0>1> 


J 


(ID 
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Notice  in  equations  (ll)  for  supersonic  speeds  that  when  the  trigono- 
metric terms  are  expanded  for  the  slow  frequency  case,  as  was  done  in 
the  previous  example,  there  appears  an  in-phase  term  of  zero  order 
in  k,  oC^aa),  and  an  out-of-phase  term  of  first  order  in  k, 

9 ■ OgL 

- ,2£_  J F1(<p)d<p.  These,  then,  are  the  principal  contributions  to  the 

lift  forces  for  the  slow-frequency  angle -of -attack  variation.  The  phase 
relationships  for  this  motion  are  shown  graphically  in  sketch  (i).  For 
this  case,  it  is  evident  that  the  total- lift  force  can  lag  behind  the 
angle  of  attack  because  of 
the  negative  out-of-phabe 
contribution, 

— f a F,(<p)  cos  k<pdq> 

AtT  A 


2V0  'o 


Application  of  indicial 
functions  to  harmonic  rotary 
oscillations.-  Finally,  con- 
sider  the  case  of  harmonic 
rotary  oscillations.  Here, 
as  previously  mentioned,  the 
normal  velocity  over  the 
wing  surface  is  composed  of 
contributions  from  both  the 
angular  velocity  and  the 
instantaneous  angle  of 
attack,  so  that  the  complete 

expression  in  equation  (3)  roust  be  employed  to  obtain  the  total  lift. 
However,  for  single-degree-of-freedom  rotary  oscillations,  a equals  9 
and  d equals  q,  so  that  in  this  case  the  separate  expressions  given 
for  the  harmonically  pitching  wing  (eqs.  (7)  and  (8))  and  the  harmoni- 
cally plunging  wing  (eq.  (ll))  can  be  combined  to  give  the  total  lift 
for  a wing  executing  harmonic  rotary  oscillations.  Then,  adding  the 
results  of  equations  (7)>  (8 ),  and  (ll),  the  in-phase  and  out-of -phase 
components  of  the  total  lift  become 


^ = Ci^U)  - k /X  F1(q>)  sin  kq>d<p  + k2/X  Fa(q>)  cos  k*Pd<P 


(12a) 


X = «B,  Mq<  1 
k~Oa>  Mq>  1 


q8/2Vc 


CLq(k)  - k £ F2(<P)  sin  WPd<p  - / FX(<P)  cos 


k<Pd<P 


(12b) 
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The  phase  relationships  for  the  rotary  oscillation  may,  of  course,  also 
be  obtained  by  directly  adding  the  results  given  in  sketches  (h)  and  (i) 
The  result  of  this  addition  is  shown  in  sketch  (j). 


It  will  be  noted  in  sketch  (j)  that  the  total-lift  force  can  either  lag 
behind  or  lead  the  angle  of  attack,  depending  on  the  relative  magnitudes 
of  the  three  terms  comprising  the  out-of-phase  lift.  The  total  lift  is 
shown  lagging  behind  the  angle  of  attack  in  sketch  (j),  which  situation, 
for  axis  positions  ahead  of  the  point  of  concentration  of  the  total  lift 
gives  rise  to  the  possibility  of  the  development  of  negatively  damped 
rotary  oscillations. 

Again,  the  complete  frequency-dependent  equations  for  the  total 
lift  of  a wing  in  supersonic  flow  due  to  the  rotary  oscillation 
(eqs.  (12),  X ■ oa)  may  be  reduced  to  first  order  in  k for  the  slow- 
frequency  case  in  the  same  manner  as  was  described  in  the  two  previous 
examples  to  give 


T " Clo(o«) 

l (13) 

Wkl  ■ CM°*>  - C Fl(<P)49 
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For  all  three  examples,  the  same  procedure  may,  of  course,  be  used 
to  obtain  the  pitching-moment  coefficient.  Only  the  pitching -moment 
equations  for  the  rotary  oscillation  are  presented  below,  since  the 
correspondence  between  the  lift  and  moment  equations  is  obvious.  For 
the  rotary  oscillation  case,  then, 

= cmo.(l)  - k fX  FS(<P)  sin  k<pd<p  + k2  fX  F4(<P)  cos  k<Pd<P 
CL  ^ o O 

\=»,  Mq< 
l=0a>  Mo> 

-^rQ  - C^Aj-k  £ F4(<p)  Sin  kq>d<p  - fx  F3(<P)  cos  kqpd<p 


where,  as  previously, 

Fa(<p)  = C^x)  - 0^(9) 

and 

F4(<P)  = Cmq(\)  - 0^(9) 

Again,  reducing  equations  (14)  for  supersonic  speeds  (x=oa)  to  first  order 
in  frequency,  there  results 


The  complete  frequency-dependent  equations  for  the  lift  and  pitching- 
moment  coefficients  for  the  rotary-oscillation  case  (eqs.  (12)  and  (l4)) 
describe  completely  the  aerodynamic  forces  and  moments  resulting  from 
the  single-degree-of-freedom  pitching  mode.  For  the  purposes  of  the 
present  discussion,  however,  it  is  sufficient  to  limit  consideration  to 
the  simpler  first  order  in  frequency  results  of  equations  (13)  and  (15). 
The  significance  of  the  effect  of  the  higher-order  terms  on  the  out-of- 
phase  pitching-moment  will  be  examined  in  a later  section  of  this  report. 

With  regard  to  the  first-order  results  for  the  supersonic  lift  and 
pitching-moment  coefficients,  it  is  instructive  to  note  that  the 
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ua  ua 

quantities  f F,(9)d<p  and  / F«((p)dq>  in  equations  (13)  and  (15) 

are  represented  geometrically  by  the  areas  of  the  shaded  portions  of 
sketch  (k) . 


Sketch  (k) 


The  manner  in  which  these  areas  are  affected  by  variations  in  Mach 
number,  aspect  ratio,  plan-form  shape,  and  thickness  will  be  used  as  a 
guide  in  later  sections  of  this  report  to  determine  the  significance  of 
these  parameters. 

In  the  foregoing  discussion,  no  mention  has  been  made  of  reducing 
the  complete  equations  for  the  lift  and  moment  coefficients  at  subsonic 
speed  (eqs.  (12)  and  (14),  X=*)  to  first  order  in  frequency  as  was  done 
for  the  equations  noted  as  applying  at  supersonic  speed.  It  is  evident 
that  if  the  same  procedure  had  been  applied  for  subsonic  speeds,  the 
area  corresponding  to  the  term  / Fl(<P)d<P  can  either  be  finite  or  can 
become  infinitely  large,  depending  on  the  manner  in  which  the  indie ial 
lift  function  Cj^^cp)  approaches  its  steady-state  asymptote  as 
In  the  latter  case,  there  exists  the  interesting  anamoly  of  an  infinite 
out-of-phase  lift  force  as  the  frequency  approaches  zero.  As  can  be 
seen  from  the  results  of  reference  4,  such  will  be  the  case  for  the  two- 
dimensional  wing.  This  result  as  the  frequency  approaches  zero  is  not 
peculiar  to  the  indicial  analysis  alone,  but  has  been  pointed  out  by  a 
number  of  authors  using  different  approaches.  As  indicated  by  Miles 
in  reference  11,  however,  the  anamolous  result  can  be  considered  to  be 
a consequence  of  assuming  a two-dimensional  flow,  and  there  is  reason 
to  believe  that  the  difficulty  as  the  frequency  approaches  zero  will  not 
exist  for  finite-span  wings. 

As  has  been  mentioned  previously,  the  use  of  the  characteristic 
areas  / a Fi^dcp  and  /°a  F3(9)dq>  will  be  shown  to  be  of  considerable 
value  in  estimating  the  damping- in-pitch  characteristics  of  wings  at 
supersonic  speeds.  For  the  two-dimensional  wing  at  subsonic  speeds,  the 
singularity  as  k-*0  prevents  the  use  of  such  a simplified  approach 
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without  further  study.  However,  rather  than  return  to  the  use  of  the 
full  frequency-dependent  equations  ((12)  and  (l4)),  the  reduction  of 
the  equations  for  the  out-of-phase  lift  and  moment  to  first  order  in 
frequency  will  be  made  in  such  a manner  as  to  preserve  the  significance 
of  these  areas.  To  accomplish  this  end,  equation  (12b)  is  reconsidered. 
It  is  evident  that  the  first  integral  in  equation  (12b)  may  be  discarded, 
since  its  contribution  to  the  out-of-phase  lift  is  at  least  of  second 
order  in  frequency.  The  second  integral  is  divided  into  two  parts: 


f F!(cp)cos  kCpdtp  = Fi(cp)cos  kcpdcp  + / Fx(9)  cos  kcpdcp 

Jo  o 

where  is  chosen  such  that  Fi(<Pi)  is  close  to  zero.  The  first 

integral  in  equation  (l6),  being  bounded,  then  causes  no  difficulty. 
Expanding  the  trigonometric  term  and  retaining  only  the  first  term  in 
the  expansion,  there  results  f ^ Fx ( <p) (d<p,  which  is  the  characteristic 

area  out  to  the  point  9i«  Now  for  large  values  of  9,  F1(9)  is  approx- 
imated in  reference  4 by 


(16) 


Fi(<P) 


U{a  + (p 


b '] 

+ i^9p; 


where  the  values  of  p,  a,  and  b are  dependent  on  Mach  number,  and  are 
given  for  M0  = 0,  0.5/  and  0.8  in  reference  4.  Inserting  this  quantity 
in  the  second  integral  in  equation  (l6),  we  have 


,•„  /■•cos  k9  ^ , r*  cos  k 

fx(«oob  vm  = » rTT'  a<l>  + ^ 1^51 


cos  k9 


d9 


(17) 


Performing  the  indicated  integrations  in  equation  (17),  there  is  obtained 
a term,  pb/a+91,  from  the  second  integral,  and  a term,  -pCi-|k[a+91  ]j> 
from  the  first  integral.  For  small  values  of  the  argument,  the  cosine 
integral  is  approximated  by  (see  ref.  12), 


Ci  ‘^k[a+91]l  « in-jrkfa-^jj 
k[a+9!]— > 0 

wnere  y is  Euler's  constant,  1.78107* 

Then,  through  the  first  order  in  frequency. 


cl 


■ 9i 


qc/2Vc 


= CLq(»)  - / F1(9)*P  + G(9x,k) 


(16a) 
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where 


G(<Pi,k) 


U Zn-j^k  [a+Cpi  ] j-  - 


vb 

a+q?! 


The  out -of -phase  moment  for  pitching  about  the  leading  edge  follows  from 
the  above  development,  with  the  added  result  of  reference  4 that  for 
large  values  of  <p, 


F3(<p)  « 


Fi(<P) 

4 


Then,  through  the  first  order. 


Cm 

qc/2v0 


(18b) 


Thus,  after  fixing  <Pi,  choosing  a (small)  value  of  k,  and  computing 

G(<Pi  ,k) , the  finite  areas  corresponding  to  the  terms  - /<^1F1(9)d9 
q?,  o 

and  - f F3(9)d9  can  be  assessed  in  the  same  manner  as  will  be  done 

Jo 

for  the  supersonic  case.  The  advantages  of  such  a procedure  will  be 
evident  later. 


Correspondence  between  indicial  lift  and  moment  analysis  and  other 
methods . - Before  proceeding  further  with  applications  of  the  indicial 
response  method,  it  is  appropriate  to  discuss  the  relationship  of  this 
approach  to  other  widely  used  methods. 

Following  the  fundamental  papers  of  Bryan  and  Routh,  which  intro- 
duced the  basic  differential  equations  of  motion  of  rigid  bodies  and 
their  stability  criteria,  the  historical  development  of  the  theory  of 
longitudinal  motions  of  an  aircraft  evolved  separately  in  two  fields  of 
research:  dynamic  stability  and  flutter.  Workers  in  the  dynamic  stabil- 
ity field  soon  found  that  the  longitudinal  oscillations  of  a rigid  air- 
craft in  flight  were  generally  of  small  reduced' frequency.  On  this 
basis,  the  constants  due  to  the  aerodynamic  properties  of  the  airframe 
which  appear  in  the  differential  equations  of  motion  were  considered  to 
be  independent  of  frequency.  As  a first  approach  to  the  problem  of 
obtaining  the  necessary  aerodynamic  coefficients  analytically,  the 
instantaneous  normal  velocity  distribution  at  the  surface  of  the  airfoil 
was  assumed  to  be  constant  with  time.  The  aerodynamic  forces  and  moments 
arising  from  the  fixed  boundary  conditions  were  then  calculated  using 
steady-flow  theory.  Later,  this  assumption  was  realized  to  be  an  over- 
simplification for  the  case  of  wing-tail  combinations  and  an  additional 
term  correct  to  the  first  order  in  frequency  was  added  which  accounted 
for  the  lag  in  the  tail  pitching  moment  caused  by  the  time  required  for 
the  vorticity  discharge  from  the  wing  to  reach  the  tail  (see  ref.  13). 
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Since  at  low  speeds  the  pitching  moment  of  the  tail  far  outweighs  all 
other  contributions,  the  results  from  steady- flow  theory  together  with 
the  term  accounting  for  the  vorticity  lag  satisfactorily  predicted  the 
dynamic  longitudinal  motions  of  wing-tail  combinations,  and  it  was  con- 
cluded that  the  major  aerodynamic  effects  had  been  accounted  for.  In 
recent  years,  however,  numerous  authors  (in  particular,  Miles;  see,  e.g., 
ref.  11)  have  pointed  out  that  the  above-mentioned  theory  overlooks 
important  contributions  to  the  aerodynamic  forces  and  moments  which, 
though  still  within  the  first  order  in  frequency  approximation,  arise 
from  time-dependent  boundary  conditions  and  must  be  calculated  from 
unsteady-flow  theory.  It  has  been  shown  by  these  authors  that  with 
proper  inclusion  in  the  equations  of  motion  of  these  coefficients,  the 
deterioration  of  damping  in  the  short-period  mode  actually  occurring  for 
aircraft  flying  at  speeds  near  the  speed  of  sound  can  be  successfully 
predicted.  The  consequences  of  the  assumptions  involved  in  the  classical 
dynamic  stability  theory  will  be  more  evident  from  a brief  review  of  the 
equation  of  motion  and  boundary  conditions  for  the  single-degree-of- 
freedom  rotary  oscillations  of  a rigid  wing  flying  at  constant  supersonic 
speed.  At  the  very  outset,  the  assumption  is  generally  made  that  the 
aerodynamic  reactions  to  the  motion  of  the  airframe  depend  only  on  the 
angular  position  and  angular  velocity  and  not  upon  angular  accelerations 
or  higher  time  derivatives . The  equation  of  motion  for  the  change  in 
pitching  moment  following  a displacement  from  an  equilibrium  position  is 
then  written  in  the  form  of  a power  series: 


I .. 

- Of  ^ 

^cm  ] /w  + 

f dCm  | dc  ( 

f dCm  ) dc  ( 

d2Cm 

{JL,  — • 

qQSc 

- da- 

Id [ac/2V0  ] 1 2V0 

Id  [qc/2V0  ]J  2V0 

_ d2a- 

J 

r £>2cm 

1 [dc/2V0]2  f 

d2Cm  1 [dS/2V 

sC  + 

1 

U2[dc/2V0] 

|J  21  t 

d2[q5/2VQ]J  21 

[AoT 

21 


(19) 


It  should  be  remembered  that  for  the  rotary-oscillation  case,  the  air- 
foil is  subjected  to  changes  in  both  angle  of  attack  a and  angular 
velocity  q,  and  that  these  motions  produce  normal  velocity  patterns  at 
the  airfoil  surface  which  are  different  in  character.  Thus,  although 
for  the  single-degree-of-freedo®  case,  a and  q are  equal,  nevertheless 
their  separate  effects  must  be  considered  and  it  is  therefore 


necessary  to  include  both 


dCm 
tqc/2V0] 


and 


dC, 


m 


d[dc/2V0J 


in  equation  (19)* 


Next,  if  it  is  assumed  that  the  moments  are  linearly  dependent  on  their 
respective  variables,  the  higher-order  terms  in  equation  (19)  may  be  dis- 
carded and  the  remaining  partial  derivatives  considered  as  constants  for 
the  given  wing.  There  remains,  therefore,  a linear  second-order  system 
with  constant  coefficients.  In  order  to  calculate  the  coefficients 
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(termed  stability  derivatives)  theoretically  it  became  necessary,  for 
lack  of  more  refined  theoretical  methods,  to  assume  that  the  instanta- 
neous normal  velocity  of  the  flow  at  the  surface  of  the  wing  was  fixed 

&Cm 

with  respect  to  time.  Thus,  the  partial  derivative  S could 

Mqc/2V0] 

be  calculated  as  the  pitching  moment  due  to  a constant  pitching  rate, 
that  is,  Cjn^Og^,  while  the  derivative  dCm/da  becomes  the  pitching 
moment  due  to  a constant  angle  of  attack,  that  is,  Cnia(aa).  As  a con- 
sequence of  fixing  the  normal  velocity  pattern  in  time,  however,  it  was 

necessary  to  assume  that  the  derivative  ~-rr — r was  zero.  There 

d[dc/2V0] 

was  therefore  no  possibility  for  this  theory  to  predict  the  occurrence 
of  dynamic  instability  for  a wing  alone,  since  the  only  damping  term 
remaining  is  Cmq(oa)>  which  is  always  stabilizing.  When  the  restriction 
of  constant  normal  velocity  with  time  is  lifted,  however,  the  assumption 
is  then  made  that  the  stability  derivatives  in  equation  (19)  may  be  calcu- 
lated separately  by  fixing  each  of  the  independent  variables  a,  a, 
and  q in  turn  with  respect  to  time.  The  derivatives  Cmq(aa) 

and  C™  (ca)  thus  remain  unchanged,  but  the  derivative  ■ — (or  Cm*) 


d[dc/2V0] 


can  now  be  included  and  calculated  as  the  pitching  moment  due  to  a con- 
stant vertical  acceleration,  aV0.  It  should  be  emphasized  that  while 
Cniq(oa)  and  Cma(oa)  may  be  calculated  from  steady- flow  theory  by  virtue 

of  the  assumed  invariance  with  time  of  the  normal  velocity  pattern, 
must  be  calculated  from  unsteady-flow  theory  since  for  constant  & 
the  angle  of  attack  varies  linearly  with  time,  as  does  the  normal  veloc- 
ity of  the  flow  at  the  surface. 


It  is  clear  that  since  the  stability  derivatives  in  equation  (19)  are 
assumed  to  be  independent  of  the  frequency,  the  result  for  the  aerody- 
namic pitching-moment  coefficient  is  thereby  limited  to  one  that  is 
correct  only  to  the  first  order  in  frequency.  For  the  single-degree-of- 
freedom  case,  then,  the  in-phase  and  out-of -phase  components  of  the 
total  aerodynamic  pitching-moment  coefficient  correct  to  the  first  order 
in  frequency,  become 


Cm 

"a" 


1 

* 


~ = Cm^(oa)  + Cm4(oa) 


(20) 


By  comparison  with  the  first  order  in  frequency  result  from  the  indicial 
response  analysis  (eq.  (15))>  it  is  evident  that  the  two  results  are 
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identical  if  the  quantity  - /°a  F3(<p)d(p  can  be  shown  to  be  equivalent 

o 

to  Cm^cTa)*  To  show  this  equivalence,  consider  a wing,  initially  in 
level  steady  flight,  which  is  suddenly  forced  down  with  constant  vertical 
acceleration  d,V0.  As  seen  in  sketch  (l),  the  angle-of-attack  variation 
in  this  case  is  a = at,  where  & is  a constant.  Then  applying  the 
counterpart  of  equation  (3)  for  the 


and  let  t be  greater  than  ta. 

Then  Sketch  (l) 


Cm(t)  = aCtKmJta)  - a / ta  F3(T)dT 

and  nondimens ionali zing,  by  replacing  t and  ta  by  c<p/2V0  and  coa/2V0, 
we  have 

Cm(<P)  = aWCmJioa)  - — /°a  F3(<p)d<P;  cp>oa  (21) 

2V0  ° 

Thus,  the  pitching  moment  proportional  to  the  constant  vertical  acceler- 
ation parameter  dc/2V0,  which  is  synonomous  with  the  definition  of  the 
stability  derivative  Cn^,  is  found  to  be  equivalent  to  the  pitching- 
moment  contribution  due  to  a for  the  first  order  in  frequency-rotary- 
oscillation  case.®  Therefore,  the  results  of  the  indicial  response 
method,  when  reduced  to  the  first  order  in  frequency  for  supersonic 
speeds,  are  identical  to  the  results  from  the  familiar  first-order  theory 
used  in  dynamic  stability  work.4 

®By  the  same  procedure,  the  stability  derivative  CL<i(aa)  can  be  shown 
to  be  equivalent  to  the  term  - C Fi(*>d<p- 
4Notice  in  equation  (21)  that  if  oa  is  replaced  by  infinity,  the  results 
apply  to  subsonic  speeds.  For  the  two-dimensional  wing,  the  analogy 

between  and  - f Fa(q>)d<p  then  gives  only  the  previously  mentioned 
singularity  at  infinity  as  k— > 0.  If  the  area  corresponding  to 
/ F3(<P)d<P  were  finite,  however,  the  analogy  would  be  equally  useful 
for  sub sonic  as  well  as  supersonic  speeds. 
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Workers  in  the  field  of  flutter,  who  were  concerned  with  frequencies 
many  times  those  encountered  in  dynamic  stability  analyses,  required 
theoretical  information  showing  the  behavior  of  the  forces  and  moments 
as  affected  by  the  frequency  of  oscillation,  and  therefore  discarded  the 
first-order  theory  for  more  precise  methods  of  analysis.  One  of  the 
most  useful  of  these  has  been  the  "oscillating  potential"  theory,  which 
is  based  on  periodic  solutions  of  the  time -dependent  linearized  equation 
of  compressible  potential  flow.  The  in-phase  and  out-of-phase  lift  and 
moments  are  thereby  determined,  generally  as  functions  of  powers  of  the 
reduced  frequency,  aspect  ratio,  Mach  number,  and  position  of  the  axis 
of  rotation.  The  application  of  this  method,  which  developed  primarily 
as  a result  of  Theodorsen's  work  for  incompressible  flow  (ref.  14)  has 
recently  produced  a number  of  useful  papers  covering  a wide  variety  of 
wings  at— supersonic  speeds  (see,  e.g.,  refs.  15  to  19) • It  has  been 
shown  by  a number  of  authors,  in  particular,  Garrick,  in  reference  7, 
that  through  the  use  of  superposition  methods  the  results  for  the  aero- 
dynamic coefficients  obtained  from  the  oscillating  potential  theory  are 
wholly  compatible  with  those  of  the  indicial  response  method  for  har- 
monic motions  (eqs.  (12)  and  (l4)). 

Thus  the  indicial  response  method  embraces  both  the  first-order 
theory  of  dynamic  stability  and  the  oscillating  potential  theory,  and, 
in  effect,  bridges  the  gap  between  the  fields  of  dynamic  stability  and 
flutter.  Furthermore,  the  indicial  response  method  overcomes  the  dis- 
advantages of  the  two  methods  described  above  since,  unlike  the  first- 
order  theory,  the  effects  of  frequency  on  the  aerodynamic  coefficients 
can  be  determined,  and,  unlike  the  oscillating  potential  theory,  the 
method  can  be  applied  easily  to  the  study  of  arbitrary  motions.  Finally, 
the  indicial  response  method  represents  a fundamental  approach  to  the 
problem  of  unsteady  flow,  and  affords  valuable  insight  into  the  physical 
nature  of  the  aerodynamic  phenomena  taking  place. 

Physical  concepts  relating  to  the  indicial  loading.-  It  has  been 
shown  that  for  even  small  frequencies,  the  pitching  moment  of  an  airfoil 
in  harmonic  rotary  motion  can  lag  behind  the  angle  of  attack  of  the  air- 
foil. The  magnitude  of  the  lag  depends  on  the  character  of  the  indicial 
response  to  a step  change  in  angle  of  attack.  It  is  therefore  of 
interest  at  this  time  to  re-examine  the  physical  nature  of  the  flow 
that  contributes  this  lag. 

Consider  first  the  lift  and  moment  at  the  instant  the  angle  of 
attack  changes,  assuming  that  previous  to  time  zero  the  wing  has  been 
flying  a level  path  at  zero  angle  of  attack.  At  t = 0,  the  wing  begins 
to  sink,  without  pitching,  with  constant  downward  velocity  aVQ  while 
maintaining  its  forward  velocity.  The  angle  of  attack  therefore  changes 
discontinuously  from  zero  to  a constant  a.  At  the  same  instant,  the 
step  change  in  angle  of  attack  causes  the  emission  of  a compression  wave 
from  each  point  on  the  lower  surface  of  the  wing  and  expansion  waves 
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from  points  on  the  upper  surface.  In  the  infinitesimal  time  during  which 
the  starting  action  occurs,  each  section  of  the  wing  experiences  the 
same  impulsive  force,  and  by  equating  the  impulse  to  the  momentum  trans- 
mitted to  the  mass  of  fluid  affected  by  the  starting  waves,  the  starting 
lift  coefficient  can  easily  be  derived  as  4a/Mo  (see  ref.  3)*  During 
the  infinitesimal  starting  time,  the  pressure  disturbances  from  the  edges 
of  the  airfoil,  propagated  at  the  speed  of  sound,  travel  an  insignificant 
distance  and  do  not  influence  the  remainder  of  the  airfoil.  The  lift 
coefficient  is  therefore  independent  of  the  wing  plan  form.  This  remark- 
ably simple  result  for  the  starting  lift  coefficient,  which  is  valid  for 
both  subsonic  and  supersonic  speeds,  is  thus  dependent  solely  on  the 
flight  Mach  number . The  starting  pitching  moment  follows  directly  from 
the  above  result,  since  by  virtue  of  the  uniformity  of  loading  the  aero- 
dynamic center  is  located  at  the  wing  centroid  of  area. 

For  values  of  time  greater  than  zero,  however,  the  situation  differs 
radically  for  the  supersonic  and  subsonic  speed  ranges.  Consider  first 
the  supersonic  case.  As  time  passes,  the  spherical  sound  waves  emitted 
at  t = 0 grow  in  size  with  radius  a0t.  The  wing,  however,  is  moving 
forward  at  a faster  rate  than  the  rate  of  growth  of  the  starting  sound 
waves  and  thus  begins  to  emerge  from  the  influence  of  these  waves.  This 


At  t = 0,  the  starting  waves  just  cover  the  wing  and  the  loading  is 
uniform  as  described  previously.  At  t = tA,  the  starting  waves  have 
grown  in  radius  and  the  wing  has  begun  to  emerge  from  their  influence. 

On  that  portion  of  the  wing  which  has  emerged,  region  (l)  in  sketch  (m), 
the  loading  has  already  reached  its  steady-state  value.  Notice  that  in 
this  region  the  characteristic  tip  Mach  cone  has  already  formed.  On  the 
portion  of  the  wing  uninfluenced  by  the  starting  waves  from  the  edges, 
region  (3)  In  sketch  (m),  the  loading  is  still  uniform  as  at  t = 0. 

In  regions  (2)  and  (4),  the  loading  Is  Influenced  by  the  starting  waves 
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from  the  leading  and  side  edges,  and  in  these  regions  is  thus  different 
from  the  loading  in  either  regions  (l)  or  (3).  As  time  increases  still 
further,  the  uniform  starting  load  quickly  disappears  as  the  sound  waves 
from  the  leading  edge  grow  in  size  and  as  the  wing  moves  forward. 

Finally,  at  time  ta,  the  envelope  of  the  starting  waves  from  the  lead- 
ing edge  is  coincident  with  the  trailing  edge  of  the  wing,  and  the 
steady-state  loading  corresponding  to  the  new  angle  of  attack  a has 
been  completely  established  over  the  wing. 

The  above  relationships  can  be  shown  more  clearly  for  the  entire 
time  interval  zero  to  ta  for  a two-dimensional  wing  by  plotting  as  a 
function  of  time  the  position  of  the  wing  leading  and  trailing  edges 
and  position  of  the  envelopes  of  the  sound  waves  which  emanate  from  the 
leading  and  trailing  edges  at  t = 0.  Such  a plot  is  shown  in  sketch  (n). 


It  is  clear  that  at  t = t*  the  regions  of  the  wing  (l),  (2), 
and  (3)  correspond  to  the  same  regions  at  t = ti  for  the  wing  shown 
in  sketch  (m) . For  t = 0 and  in  region  (3)  the  loading  is  uniform  and 
is  given  by  4a/Mo.  For  t>ta  and  in  region  (l)  the  wing  has  out- 
stripped the  starting  waves  from  the  leading  edge  and  has  attained  its 
steady-state  loading.  For  t<ta  the  chordwise  loading  is  composed  of 
combinations  of  the  loading  in  each  of  the  three  regions  as  shown  in 
sketch  (o) .5 


5The  reader  will  note  the  similarity  between  sketch  (n)  and  sketches 
depicting  the  boundary  conditions  for  three-dimensional  wings  in 
steady  supersonic  flow.  Many  researchers  have  pointed  out  the  analogy 
and  it  has  been  used  to  calculate  the  pressure  over  a wing  impulsively 
starting  from  rest  (refs.  3 and  4). 
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It  may  be  easily  verified  that  ta,  the  time  required  for  the  wing  to 
attain  its  steady-state  loading,  is  given  by  the  following  relationship: 


cqMq 

V0[Mo-csc  v] 


(22) 


Notice  that  for  wings  having  straight  or  swept-forward  trailing  edges 
and  straight  or  swept-back  leading  edges,  equation  (22)  reduces  to 


cqMq 
V0 [Mo-1] 


(23) 


In  terms  of  the  number  of  half  M.A.C.  lengths  of  travel,  it  thus  appears 
that  for  wings  having  straight  or  swept-forward  trailing  edges  and 
straight  or  swept-back  leading  edges, 


_ 2vo  £ _ 2c0  Mq 

c a c [Mq-1 ] 


(24) 


whereas  for  wings  having  swept-back  trailing  edges  and  straight  or 
swept-back  leading  edges. 


°a  = 


oa  = 


2c0  Mq  2Cq 

— =—  p- r,  Mq  >CSC  V + — -=■  COS  V 

c [Mq-csc  v]  ^ b 


2Mof 


Mq 


2cr 


-—5  [Mo2-l]  \,  Mq^csc  v + — ^ cos  V 
4c  I b 


(25) 


Mq2  - 1 


The  second  of  equations  (25)  applies  to  that  range  of  Mach  numbers  for 
which  the  trace  of  the  starting  sound  wave  from  the  apex  is  not  tangent 
to  the  trailing  edge  at  t = ta. 


Now  consider  the  subsonic  case.  Here  the  situation  is  more  compli- 
cated in  that  since  the  starting  sound  waves  travel  faster  than  the 
wing,  the  wing  never  escapes  their  influence.  Furthermore,  the  vortic- 
ity  shed  by  the  wing  at  t = 0 can  also  influence  the  loading  on  the 
wing  since  the  disturbances  created  by  the  shed  vorticity  travel  forward 
at  a faster  rate  than  the  wing.  For  these  reasons,  the  indicial  loading 
at  subsonic  speeds  approaches  its  steady-state  distribution  asymptoti- 
cally with  time.  The  situation  for  subsonic  speeds  will  be  more  clearly 
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evident  from  examination  of  sketch  (q),  which  shows  the  relationship  of 
the  traces  versus  time  of  the  leading  ajid  trailing  edges  of  a two- 
dimensional  wing  flying  at  a subsonic  Mach  number  to  the  traces  of  the 
envelopes  of  the  starting  sound  waves  - 


Sketch  (q) 

Notice  in  sketch  (q)  that  the  starting 
sound  waves  intersect  the  edges  of  the  -wing 
and  that  each  intersection  causes  a new 
sound  wave  to  be  emitted,  which  in  turn  will 
intersect  an  edge.  Furthermore,  notice  that 
the  vorticity  shed  from  the  trailing  edge 
at  time  zero  can  influence  that  portion  of* 
the  wing  behind  the  sound  wave  trace 
labeled  a-b.  The  influence  of  each  succes- 
sive sound  wave  reflection,  however , is 
weaker  than  the  last,  and  as  the  wing  moves 
away  from  the  starting  vortices  their 
influence  diminishes,  so  that  at  t = » the 
loading  on  the  wing  attains  its  steady-state 
distribution.  The  variation  of  the  chord- 
wise  loading  with  time  for  the  two- 
dimensional  wing  flying  at  a subsonic  Mach 
number  is  shown  in  sketch  (r).  Notice  "that 
for  t >0,  the  chordwise  loading  is  markedly 
different  from  the  loading  at  supersonic 
speeds  (sketch  (o) ) . However,  for  t * 0 


Sketch  (r) 
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and  in  the  region  corresponding  to  region  (3)  of  sketch  (q),  the  loading 
is  uniform  and  equal  to  4a/Mo,  as  in  the  supersonic  case. 


Damping  in  Pitch  of  Low-Aspect-Ratio  Wings 


Previously  (eq.  (20)),  it  was  shown  that  for  single-degree-of- 
freedom,  slow-frequency,  pitching  oscillations  of  a wing,  the  principal 
parameter  contributing  to  the  damping  of  the  motion  is  the  damping  coef- 
ficient Cmq  + Cju^.  This  result,  however,  is  not  directly  applicable  to 
the  case  of  an  aircraft  in  flight,  since  generally  additional  damping  is 
provided  by  virtue  of  the  fact  that  the  aircraft  experiences  harmonic 
vertical  translatory  oscillations  as  well  as  the  rotary  oscillations. 

For  the  slow-frequency  case,  however,  the  effect  of  the  second  degree  of 
freedom  on  the  damping  can  be  determined  by  use  of  the  relationships 
given  in  reference  10.  It  is  further  shown  in  reference  10  that  although 
the  effect  of  the  translatory  oscillation  is  usually  to  increase  the 
total  damping,  nevertheless,  the  parameter  of  primary  importance  in 
determining  the  magnitude  and  duration  of  the  motion  remains  the  damping- 
in-pitch  coefficient  + Cnu.  The  remainder  of  this  section  is 

therefore  concerned  witha  study  of  the  effect  on  this  parameter  of 
certain  important  variables.  In  particular,  the  effect  of  the  position 
of  the  center  of  gravity,  and  the  effects  of  Mach  number,  aspect  ratio, 
plan-form  shape,  frequency,  and  thickness  will  be  examined,  principally 
by  inspection  of  the  indicial  lift  and  moment  responses  to  a change  in 
angle  of  attack. 

Effect  of  static  margin.-  From  the  previous  discussion,  it  will  be 
remembered  that  at  supersonic  speeds  the  stability  derivatives  Cp^(oa) 
and  Cnu(oa)  were  shown  to  be  equivalent  to  the  indicial  lift  and  moment 

CJgL  0 

expressions,  - f F1(qp)d9  and  - / F3(cp)dcp.  Furthermore,  it  was  shown 

that  these  were  the  contributions  which  could  cause  the  total  lift  and 
moment  during  the  short-period  oscillation  to  lag  behind  the  angle  of 
attack.  Hence,  by  inspection  of  equation  (15),  it  is  evident  that  since 
Cmq(°a)  Is  always  stabilizing,  when  Cma(0a)is  negative  (corresponding 
to  a statically  stable  condition)  the  possibility  of  dynamic  instability 
in  the  form  of  negatively  damped  rotary  oscillations  arises  when 
Cm&(aa)>  the  shaded  area  in  sketch  (k-ii),  is  larger  than  (^(oa,). 

Now  since  the  normal  velocity  at  the  surface  of  the  wing  due  to  the 
instantaneous  angle  of  attack  is  constant  over  the  wing,  the  lift 
derivative  Cl^  which  arises  from  this  boundary  condition  is  independ- 
ent of  the  axis  position  and  will  therefore  vary  linearly  with 
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axis  position.  This  variation  is  illustrated  in  sketch  (s),  where  Axo 
represents  the  distance  of  the  center  of  gravity  from  the  aerodynamic 
center,  measured  positive  forward  of  the  aerodynamic  center. 


The  parameter  CLq,  on  the  other  hand,  is  a direct  function  of  the 
axis  position,  since  it  arises  from  a normal  velocity  distribution  that 
varies  directly  as  the  distance  from  the  axis.  The  moment  coeffi- 
cient Cniq  will  therefore  vary  as  the  second  power  of  the  axis  posi- 
tion, and  describes  the  parabolic  shape  shown  also  in  sketch  (s).  It 
is  evident  from  sketch  (s)  that  the  sum  of  Cmn  and  will  be  a 

minimum  at  some  value  of  the  static  margin,  ana  that  the  sign 
of  Cmq  + Cm^  at  that  point  determines  whether  or  not  a region  of  axis 
positions  will  exist  over  which  the  wing  can  experience  negatively 
damped  rotary  oscillations.  These  qualitative  statements  may  be  written 
explicitly  by  considering  the  equation  for  the  damping  in  pitch  about  an 
arbitrary  axis. 


where  again  Ax<>  refers  to  the  distance  of  the  center  of  gravity  frc-* 
the  aerodynamic  center,  and  the  subscripted  terms  are  referred  to  an 
axis  through  the  aerodynamic  center. 
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Taking  the  derivative  of  equation  (26)  with  respect  to  AXq/c,  and 
setting  the  result  equal  to  zero,  there  is  obtained  the  axis  location 
at  which'  the  damping  in  pitch  is  a minimum 

•\ 


d [Cmq+Omd] 
d[Ax0/c] 


0 


/ 


(27) 


When  equation  (27)  is  inserted  into  equation  (26),  the  minimum  value  of 
the  damping  in  pitch  is  given  as 


Cmq+Cma 


Jmin. 


and  hence,  a region  of  instability  will  exist  if 


1 

SCLo, 


{ 


6 


(28) 


(29) 


If  equation  (29)  is  greater  them  zero,  the  boundaries  of  the  region  of 
axis  positions  over  which  instability  is  possible  are  of  course  given  by 
setting  equation  (26)  equal  to  zero  and  solving  for  AXq/c. 


[.cLq)  +cLi' 


'}  + { 


[cmqJ0+ 


(30) 


Notice  in  both  sketch  (s)  and  equation  (30)  that  for  a given  Mach  number 
there  will  be  two  axis  positions  at  which  the  deunping  in  pitch  vanishes. 
Then  if  the  above  procedure  is  carried  out  for  a series  of  Mach  numbers, 
one  may  trace  out  a curve  as  shown  in  sketch  (t)  which  forms  the  locus 
of  Mach  numbers  and  axis  positions  at  which  the  damping  in  pitch  is  zero. 
This  locus  thus  delineates  the  regions  of  Mach  number  and  axis  position 
for  which  dynamic  instability  is  and  is  not  theoretically  possible. 

Such  loci,  covering  a wide  variety  of  wings  at  supersonic  speeds,  have 
been  presented  by  a number  of  authors.  Watkins,  for  example,  presents 
supersonic  boundary  curves  for  rectangular  and  triangular  wings  in  refer- 
ences 16  and  17.  At  subsonic  speeds.  Miles*  reduction  of  Possio's 
development  to  first  order  in  frequency  (ref.  11)  can  be  used  to  form  a 
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stability  boundary  curve  for  the  two- 
dimensional  wing  for  a given  (small) 
reduced  frequency.  Such  a curve  is  pre- 
sented for  the  entire  Mach  number  range 
in  figure  1,  where,  here  + xQ  is  the 
distance  of  the  axis  of  rotation  behind 
the  leading  edge  and  k,  the  reduced  fre- 
quency, is  0.011  for  subsonic  speeds  and 
approaches  zero  for  supersonic  speeds. 

Notice  in  figure  1 that  at  both  subsonic 
and  supersonic  speeds,  the  range  of  Mach 
numbers  over  which  dynamic  instability  is 
possible  is  largest  for  center-of-gravity 
positions  forward  of  the  aerodynamic 
center.  Further,  the  largest  range  of 
axis  positions  over  which  dynamic  insta- 
bility is  possible  occurs  near  Mq  = 1. 

Both  of  these  characteristics  are  true  as  well  for  three-dimensional 
wings  at  supersonic  speeds  (see,  e.g.,  refs.  16  to  19). 


Effect  of  Mach  number.-  Next  consider  the  effect  of  Mach  number  on 
the  damping  in  pitch  of  a two-dimensional  wing  with  axis  at  the  leading 
edge.  The  variation  with  Mach  number  of  the  indie ial  pitching-moment 
response  to  a change  in  a will  first  be  examined,  using  the  information 
given  in  the  previous  sections  and  the  indicial  curves  given  in  refer- 
ence 4.  At  supersonic  speeds,  the  manner  in  which  - fo°aF3(<p)d<pf  the 


area  corresponding  to  Cm£,  is  affected  can  then  be  assessed  and  com- 
pared with  Cjnn(ca)*  At  subsonic  speeds,  use  is  made  of  equation  (l8b). 
It  is  evident  in  equation  (l8b)  that  by  fixing  k and  choosing  <Pi> 

such  that  the  quantity  i G(cp1,k)  is  the  sane  at  each  Mach  number,  one 

4 flp. 

is  free  to  compare  finite  areas  - / Fa(<p)dq>  on  an  equivalent  basis. 


As  has  been  mentioned  previously,  the  starting  lift,  at  any  Mach 
number,  is  4a/Mo  and  Is  concentrated  at  the  midchord.  At  Mo  = 0, 
therefore,  there  is  an  initial  infinite  pulse  in  the  pitching  moment 
about  an  axis  coincident  with  the  leading  edge  after  which  the  indicial 
curve  drops  to  n/4  and  begins  to  grow  asympotically  toward  its  steady- 
state  value  jt/2.  At  low  subsonic  Mach  numbers  other  than  zero,  the 
initial  pitching  moment  is  less  than  infinite  but  very  large,  being 
2/Mo*  and  then  falls  before  growing  toward  the  steady  value  n/2(3 . As 
the  Mach  number  increases  toward  1.0 , the  starting  pitching  moment 
falls  while  the  asymptotic  value  grows,  until  at  Mo  ■ 1.0  the  indicial 
curve  becomes  unbounded  in  asymptotic  moment. 
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As  seen  in  sketch  (u) , the  effect  of  increasing  the  Mach  number  at 
subsonic  speeds  is  therefore  to  increase  rapidly  the  area  corresponding 

to  the  destabilizing  moment  contribution,  - Fg(<P)d<p. 


Sketch  (u) 

In  sketch  (u),  k was  chosen  to  be  0.011  and  the  values  of  9i  were 
picked  such  that  - i GC^jk)  was  +4.88  at  each  Mach  number.  In  the 
following  discussion,  the  damping  moment  - J*  Fa(<p)d<p  + 4.88  will  be 
referred  to  as  Cm^  for  convenience.  ° 


At  supersonic  speeds,  the  initial  value  of  the  pitching  moment  2/Mq 
continues  to  drop  with  increasing  Mach  number,  but  here  the  steady-state 
pitching  moment  also  begins  to  fall  and  at  a faster  rate  than  the  start- 
ing moment,  being  2/p.  Even  more  important,  as  the  Mach  number  increases, 
the  number  of  half-chord  lengths  traveled  to  reach  steady  state  decreases 
rapidly,  being  22,  for  example,  at  Mq  = 1.1,  as  compared  to  4 at  Mq  = 2. 

As  seen  in  sketch  (v),  the  area  representing  Cm^  therefore  shrinks 

rapidly  with  increasing  supersonic  Mach  number  and  becomes  relatively 
unimportant  at  Mach  numbers  greater  than  2.  The  trend  of  with 

Mach  number  through  the  range  0<Mq<2  is  more  clearly  evident  in 
sketch  (w) . It  is  seen  that  is  positive,  or  destabilizing, 

throughout  the  Mach  number  range  and  that  its  effect  is  most  important 

at  Mach  numbers  near  1.0.  Also  shown  plotted  for  comparison  in  sketch  (w) 

is  the  variation  of  Cmq  with  Mach  number.  When  the  parameters  are 
added,  it  is  evident  that  the  damping  moment  Cm^  + for  the  two- 

dimensional  wing  with  axis  at  the  leading  edge  is  destabilizing  in  the 
Mach  number  range  0<Mq<1.4i4. 
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Effect  of  aspect 
ratio.-  To  illustrate  the 
effect  of  aspect  ratio, 
it  is  convenient  to 
compare  the  supersonic 
damping-in-pitch  char- 
acteristics of  a group  of  0 

triangular  wings  having 
subsonic  leading  edges. 

The  wings  are  of  equal 
area  and  differ  only  in 
aspect  ratio.  As  was 
done  previously,  the 
inaicial  lift  responses  to 

to  a change  in  angle  of 
attack  will  first  be 

examined.  The  effect  of  _ _ 

30 

aspect  ratio  on  the  char- 
acteristic  area  repre- 
senting can  then  be 

assessed. 

As  has  been  mentioned  previously,  the  starting  lift  coefficient 
after  a step  change  in  a is  independent  of  aspect  ratio  and  is  there- 
fore equal  to  Ua/Mo  for  each  wing.  The  parameter  aa,  the  number  of 
half  M.A.C.  lengths  required  to  reach  steady  state,  is  also  the  same  for 
each  wing,  being  a function  only  of  Mach  number.  The  steady-state  lift 
coefficient,  on  the  other  hand,  is  a function  of  aspect  ratio  and 


Sketch  (w) 
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decreases  as  the  aspect  ratio  is  reduced.  Thus,  as  shown  schematically 
in  sketch  (x),  as  the  aspect  ratio  becomes  smaller,  the  characteristic 


Sketch  (x) 


area  representing  decreases  rapidly.  For  the  wing  of  smallest 

aspect  ratio,  Cl£  may  be  positive  since  the  area  below  the  steady- 
state  lift  coefficient  is  more  than  compensated  for  by  the  area  above. 

It  is  evident,  therefore,  that  a reduction  in  aspect  ratio  has  a highly 
stabilizing  effect  on  the  damping  in  pitch,  since  for  positive  values 
of  the  static  margin  the  development  of  a destabilizing  damping  moment 
is  possible  only  when  is  negative.  This  result  is  shown  in 

sketch  (y)  where,  for  an  axis  of  rotation  located  at  0.25  c and  Mq=1.2, 
the  damping  parameters  are  presented  as  functions  of  aspect  ratio. 

Since  for  triangular  wings  the  lift  due  to  4 is  concentrated  at  % 5, 


is  equal  to  - ^ C]^. 


The  variation  of  Cjd»  with  aspect  ratio 


g '-Ha'  ““a 

shown  in  sketch  (y)  then  follows  directly  from  the  trend  of  Cj^  shown 


in  sketch  (x).  Also  plotted  in  sketch  (y)  is  the  variation  of  G 


mq 


eTheoretical  indicial  curves  have  not  yet  been  calculated  for  the 
triangular  wing  with  subsonic  leading  edges.  The  curves  drawn  in 
sketch  (x)  are  estimates  of  the  true  shapes,  and  are  intended  only 
to  indicate  the  trend  of  the  characteristic  area  with  aspect  ratio. 
The  exact  variation  of  Cjj;  with  aspect  ratio  can  be  computed  from 
the  results  of  reference  20. 


CONFIDENTIAL 


NACA  RM  A52L04a 


CONFIDENTIAL 


35 


Sketch  (y) 


with  aspect  ratio  (ref.  20).  It  is  apparent  that  although  Cm^  becomes 
more  stabilizing  with  increasing  aspect  ratio,  the  destabilizing  effect 
of  predominates,  and  the  trend  of  the  net  damping  moment  is  seen 

to  become  highly  destabilizing  as  the  aspect  ratio  is  increased. 

By  the  same  reasoning,  the  variation  with  aspect  ratio  of  the 
damping  moments  of  other  types  of  wings  can  be  shown  to  be  similar  (see, 
e.g.,  refs.  11,  21,  22,  and  23).  A notable  exception,  however,  is  the 
triangular  wing  with  supersonic  leading  edges,  whose  damping  in  pitch 
has  been  shown  to  be  independent  of  aspect  ratio  (see  refs.  18  and  19)* 
This  characteristic  may  be  anticipated  from  a study  of  the  indicial 
response  curves,  since  not  only  are  the  initial  pitching  moment  Cm^o) 
and  the  half  M.A.C.  lengths  traveled  to  reach  steady  state  (oa)  inde- 
pendent of  aspect  ratio,  but,  unlike  the  sub sonic -edged  triangular  «d.ng, 
the  steady-state  pitching  moment  Cma(oa)  is  also  independent  of  aspect 

ratio,  being  - — — £ L Inspection  of  the  results  of  reference  4 then 

P L J 

reveals  that  the  indicial  variation  Cma(  9)  between  zero  and  oa  and 
the  steady-state  parameter  Cm  (ca)  are  likewise  independent  of  aspect 
ratio. 

Effect  of  plan-form  shape.-  Next,  consider  the  effect  of  plan-form 
shape  on  the  supersonic  damping  in  pitch  of  a group  of  wings  having  the 
same  aspect  ratio.  In  order  to  make  use  of  readily  available  theoretical 
results  and  yet  compare  these  results  with  the  experimental  data  to  be 
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given  later,  three  wings  of  aspect  ratio  3 are  chosen,  having  the  trian- 
gular, swept,  and  rectangular  plan  forms  shown  in  sketch  (z). 


Sketch  ( z) 

As  has  been  shown  in  the  section  entitled  "Effect  of  static  margin," 
the  tendencies  of  the  wings  toward  dynamic  instability  can  be  compared 
comprehensively  by  plotting  their  stability  boundaries.  For  this  com- 
parison, then,  use  is  made  of  equation  ( 30) . The  stability  derivatives 
which  appear  in  equation  (30)  were  computed  from  the  theoretical  results 
of  references  11,  20,  22,  and  23.  Results  of  these  calculations  are 
shown  in  figure  2,  wherein  the  stability  boundaries  for  the  three  wings 
are  shown  as  a function  of  axis  position  and  Mach  number.  (Note  that 
the  axis  position  for  each  wing  is  measured  as  the  distance  from  the 
leading  edge  of  the  M.A.C.  of  the  wing,  and  that  the  dimensions  are  non- 
dimensionalized  on  an  equivalent  basis  by  referring  them  to  the  M.A.C. 
of  the  triangular  wing.)  It  is  clear  from  inspection  of  figure  2 that 
at  any  Mach  number  the  triangular  wing  has  the  smallest  range  of  axis 
positions  over  which  dynamic  instability  is  possible  and  the  rectangular 
wing,  the  largest. 

The  differences  in  the  damping  characteristics  of  the  triangular 
and  rectangular  wings  will  be  more  clearly  understood  by  a qualitative 
study  of  their  indicial  lift  responses  for  a Mach  number  of  1.2,  and  an 
examination  of  the  distribution  of  loading  due  to  d for  the  two  wings. 
Consider  first  the  indicial  lift  responses. 

Again,  the  starting  lift  coefficient  is  independent  of  plan-form 
shape  and  is  4a/Mo  for  each  wing.  For  the  rectangular  wing,  the 
lift  drops  abruptly  after  time  zero  due  to  the  loss  in  lift  in  the 
regions  of  the  wing  influenced  by  the  formation  of  the  tip  Mach  cones 
and  the  starting  waves  from  the  side  edges  (see  sketch  (m)).  Then  as 
the  wing  begins  to  emerge  from  the  influence  of  the  starting  waves,  the 
lift  begins  to  recover,  rises  toward  its  steady-state  value  (given 
by  eq.  6.3-2  of  ref.  11),  and  attains  this  value  after  12  half-chord 
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lengths  of  travel  (eq.  (24)).  The  variation  is  shown  in  reference  4 and 
is  redrawn  in  sketch  (aa).  As  mentioned  previously,  theoretical  indicial 


Sketch  (aa) 

lift  results  have  not  yet  been  developed  for  triangular  wings  having 
subsonic  leading  edges.  However,  the  variation  shown  in  sketch  (aa)  is 
considered  to  be  a reasonable  estimate  of  the  true  shape,  being  based  on 
knowledge  of  the  steady-state  lift  (ref.  20),  the  time  to  reach  steady 
state  (eq.  (24)),  and  the  assumption  that  the  shape  of  the  variation 
would  be  similar  to  that  of  the  wide  triangular  wing  (ref.  4),  The 
curve  was  adjusted  within  the  known  end  points  until  the  area  correspond- 
ing to  agreed  with  that  given  for  this  parameter  in  reference  20. 

It  is  evident  from  examination  of  sketch  (aa)  that  because  of  the  initial 
loss  in  lift  and  the  larger  steady-state  lift  for  the  rectangular  wing, 
for  this  wing  is  significantly  more  negative  than  that  for  the  tri- 
singular  wing.  Next,  it  is  shown  in  references  11  and  21  that  with  the 
exception  of  regions  influenced  by  tip  Mach  cones,  the  loading  due  to  d 
( for  9 > oa)  for  wings  having  swept-back  leading  edges  increases  linearly 
from  zero  along  rays  from  the  apex;  whereas  for  rectangular  wings  the 
loading  due  to  A increases  linearly  from  zero  along  chord  lines.  These 
characteristics  place  the  center  of  loading  due  to  a at  — c0  for  the 

2 **  v 

triangular  wing  and  approximately  — cQ  for  the  rectangular  wing.' 

Then  for  sin  axis  of  rotation  passing  through  the  aerodynamic  centers  of 
the  wings,  the  moment  arm  for  the  lift  due  to  a for  the  triangular 


7 Due  to  the  influence  of  the  tip  Mach  cones,  the  center  of  loading  due 
to  a is  shifted  forward  somewhat  from  the  position  it  has  for  the 
two-dimensional  wing.  Calculations  for  the  A = 3 rectangular  wing 
at  Mq  = 1.2  show  that  the  center  of  loading  is  at  0.605  c0. 
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wing  is 


whereas  for  the  rectangular  wing  it  is  approxi- 


mately 


Thus,  not  only  is  the  negative  out-of-phase  lift  contri- 


bution Cl^  for  the  rectangular  wing  significantly  larger  than  that  for 

the  triangular  wing,  but  the  destabilizing  damping  moment  - £ C^i  is 

c 

larger  yet,  due  to  the  larger  moment  arm.  Calculations  for  the  steady 
pitching  parameter  Cmq  for  an  axis  through  the  aerodynamic  center 
(refs.  11  and  20)  then  reveal  that  Cm^  for  the  triangular  wing  is  more 

negative  than  for  the  rectangular  wing.  The  net  result  is  therefore  a 

considerably  larger  damping  moment  for  the  triangular  wing  than  for  the 
rectangular  wing.  The  result  of  this  comparison,  however,  should  not  be 
interpreted  as  a recommendation  that  the  triangular  rather  than  the  rec- 
tangular wing  be  used  on  aircraft  from  a dynamic  stability  standpoint. 

To  obtain  adequate  static  stability,  the  rectangular  wing  would  generally 
be  employed  in  combination  with  a tail  surface,  whereas  the  triangular 
plan  form  may  be  sufficiently  airworthy  without  the  use  of  a tail.  The 

addition  of  a tail  surface  in  effect  reduces  the  aspect  ratio  of  the 

rectangular  wing,  which  reduction,  as  noted  previously,  has  a highly 
stabilizing  effect  on  the  damping  in  pitch.  The  tailless  triangular  wing 
may  therefore  experience  more  difficulty  at  Mach  numbers  near  1.0  than  a 
rectangular  wing-tail  combination. 


Effect  of  frequency.-  The  previous  discussion  has  been  restricted 
to  the  analysis  of  a harmonic  motion  that  is  of  vanishingly  small  fre- 
quency. This  limitation  arose  as  a consequence  of  discarding  all  but 
first  order  in  frequency  terms  in  the  expansions  of  equations  (12) 
and  (l4).  The  question  arises:  when  the  frequency  can  no  longer  be  con- 
sidered small,  what  effect  has  the  frequency  on  the  damping  in  pitch? 

Previously,  the  trigonometric  terms  in  equation  ( 14)  were  expanded 
and,  assuming  k to  be  very  small,  terms  of  order  k2  and  higher  were 
eliminated.  It  was  then  found  that  the  loss  in  damping  from  that  pro- 
vided by  the  steady  damping  parameter  Cmq(aa)  was  associated  with  the 
destabilizing  moment  contribution  corresponding  to  the  term 

- J0  F3(9)d9.  Nov,  however,  we  discard  the  restriction  of  small  k 
and  perform  graphically  the  integrations  evident  in  equation  (14)  for 
supersonic  speeds  for  several  values  of  k.  The  procedure  is  indicated 
in  sketch  (bb).  It  is  apparent  from  sketch  (bb)  that  the  effect  of 
increasing  k is  to  reduce  the  area  corresponding  to  the  destabilizing 

moment  contribution  - f Fg(cp)  cos  kcpdcp.  There  appears  another 


®Again,  due  to  the  influence  of  tip  Mach  cones,  the  aerodynamic  center  is 
shifted  forward  from  c0  to  0.443  c0.  The  moment  arm  is  therefore 
0«l63  c0. 
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Sketch  (bb) 

destabilizing  contribution,  -k  £°a  F4(<p)  sin  k<pd<P,  but  quite  evidently 

its  effect  is  small  compared  to  the  reduction  in  the  term 

- /°a  Fa(<p)  cos  k<pdcp.  Notice  further  in  sketch  (bb)  that  the  effect  of 
o 

increasing  k becomes  of  marked  importance  when  the  half -period  of 
oscillation  is  the  same  order  of  magnitude  as  the  time  for  the  indicial 
response  to  reach  steady  state.  As  shown  in  sketch  (bb)  for  the  fre- 
quency k2,  the  area then  begins  to  subtract  from  (j^)>  so  that  the 
destabilizing  contribution  - f°a  F_(9)  cos  k9d<p  can  be  very  small.  We 
may  therefore  expect  that  increasing  the  frequency  of  oscillation  has  a 
stabilizing  effect  on  the  damping  in  pitch.  This  conjecture  is  sub- 
stantiated in  figures  3 and  4,  where  the  supersonic  stability  boundary 
curves  for  aspect  ratio  4 triangular  and  rectangular  wings  are  shown 
plotted  for  various  reduced  frequencies.  These  curves  were  obtained 
from  calculations  based  on  the  results  of  references  16  and  17*  Notice 
that  for  both  wings  the  region  of  possible  instability  is  diminished  as 
the  frequency  is  increased. 

From  the  results  of  the  analysis  for  supersonic  speeds,  we  may 
further  expect  that  the  stabilizing  effect  of  increasing  the  frequency 
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will  be  of  even  more  importance  at  subsonic  speeds,  for  here  the  indicial 
variation  F3(<P)  dies  out  at  infinity.  The  half-period  of  oscillation 

is  therefore  always  smaller  than  the  time  to  reach  steady  state. 

00 

The  situation  for  the  term  - J F3(T)  cos  kTd<P  is  shown  in  sketch  (cc). 


Sketch  (cc) 

It  is  evident  in  sketch  (cc)  that  the  destabilizing  moment 

- J F3(<P)  cos  kTd<P  diminishes  rapidly  as  the  frequency  is  increased. 
The  effect  of  this  reduction  on  the  damping  in  pitch  can  be  illustrated 
by  plotting  the  subsonic  damping -moment  coefficient  against  reduced 
frequency  for  the  two-dimensional  wing  (with  axis  at  the  leading  edge) 
for  Mach  numbers  0,  0.5>  and  0.8.  The  results,  which  were  obtained 
from  reference  4,  are  shown  in  figure  5.  It  is  seen  in  figure  5 that 
the  large  destabilizing  effects  of  the  moment  contributions 

- / Fa(<P)  cos  k<Pdcp  and  -k  fQ  F4(T)  sin  k<PdT  are  confined  to  a rela- 
tively narrow  range  of  reduced  frequencies.  Notice  further  in  figure  5 
that  the  range  of  frequencies  for  which  instability  is  possible  is  small 
at  Mq  = 0 (0<k<0.04)  and  grows  with  increasing  Mach  number.  This  is 
believed  to  be  the  primary  reason  why  unsteady  lift  effects  were  found 
to  be  unimportant  at  low  speedB  but  are  of  great  importance  at  speeds 
near  the  sonic  speed. 
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Effect  of  thickness.-  In  this  section,  the  concepts  regarding  the 
indicial  functions  discussed  in  previous  paragraphs  will  be  used  to 
examine  the  effect  of  profile  thickness  on  the  damping  in  pitch  of  a 
two-dimensional  wing  flying  at  a supersonic  Mach  number. 

It  will  be  recalled  that  for  low  frequencies  the  damping  in  pitch 
may  be  considered  to  be  the  sum  of  the  damping  moment  due  to  steady 
pitching,  Cm^(aa),  and  the  area  representing  on  the  plot  of  the 

pitching-moment  response  to  a step  change  in  angle  of  attack 
(sketch  (k-ii)).  The  airfoil  thickness,  of  course,  influences  both 
these  parameters.  For  the  purpose  of  this  analysis,  it  will  be  con- 
sidered sufficiently  accurate  to  study  the  effects  of  thickness  only  on 
the  steady  parameter  Cm^Oa)  and  on  the  end  values  Cji^o),  C^Ca) 
and  oa  of  the  indicial  response  curve.  It  will  then  be  assumed  that 
shape  of  the  indicial  response  curve  is  not  appreciably  altered  by  thick- 
ness effects,  whereupon  the  major  effect  of  thickness  on  Cm^  can  be 
assessed  by  adjusting  the  response  curve  given  by  the  linearized  theory 
to  fit  the  corrected  values  of  C^o),  C^Oa)  and  oa» 

It  has  been  shown  by  Busemann  (see,  e.g.,  ref.  2k)  that  the  loading 
at  an  element  of  a lifting  airfoil  may  be  represented  by 

££  = C^Ou-c^a]  + C2  jou2-oz2+2a  [ou+cj]  | (31) 

where 

Cl  =7^ 

_ [7+1]Mq4  - MMq2-1] 

2[Mo2-1]2 

and  ou  and  are  the  local  slopes  of  the  upper  and  lower  surfaces, 
respectively,  measured  with  respect  to  the  chord  line.  For  the  special 
case  of  an  airfoil  that  is  symmetrical  about  the  chord  line, 
ou  equals  Oj  and  equation  (31)  reduces  to  the  following: 

= 2aC1  + kaCg<Ju  (32) 

<Jo 

For  the  sake  of  simplicity,  in  what  follows  equation  (32)  will  be  used 
for  the  local  loading  rather  than  equation  (31) • The  results  will  there- 
fore be  valid  only  for  airfoils  that  are  symmetrical  about  the  chord  line. 
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The  lift  on  the  airfoil  due  to  angle  of  attack  is  obtained  simply  by 
integration  of  equation  (32)  across  the  chord. 


Cl/tfa) 

^La(°a) 


-L  fc°SEdx 
Co  o 

Sltsl  . 2C!  + 

a 


4C2  fcc 

— JL  cudx 


(33) 


Similarly,  the  moment  due  to  angle  of  attack,  measured  about  the  leading 
edge,  is 


Cm’  <°a)  8 * A I?  x d* 


cQ2  'o 


cma' 


(°a) 


Cm'  ( °a) 


4C. 


- - ci  - TT  4 


Co 


oux  dx 


(3M 


When  both  the  leading  and  trailing  edges  are  on  the  chord  line,  equa- 
tions (33)  and  (3*0  reduce  to  the  following: 


cLo,(  °a)  ■ 2Ci 


Cma’ 


(oa)  = - Cx  + 

co 


(35) 


where  SI  is  the  enclosed  area  of  the  airfoil  section. 

The  local  loading  given  by  equation  (32)  is  still  applicable  for 
the  steady  pitching  case  with  this  added  provision:  that  now  a repre- 
sents the  local  angle  of  attack.  For  an  airfoil  pitching  about  its 
leading  edge,  the  variation  of  a with  chordwise  distance  is  given 
by  a*  qx/v0,  where  q is  the  pitching  velocity.  Then  again,  for  the 
lift  due  to  pitching, 

♦ »• 

CL'(oa)  S^T  f ° dx 

co  o 


cLq’(°a)  * ~ = 2Ci  + ^|  /C°  ^x  dx 
^ V0/2V0  cQ*  o 


(36) 
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Likewise,  the  pitching  moment  due  to  pitching  about  the  leading  edge  is 
Cm'Ua)  = “ -ir  /C°  dx 


Co 


*0 


<V<«a)  - 1 C,  - ®2S 

<lc0/2V0  3 cca  0 


(37) 


The  primes  on  the  parameters  indicate  that  the  moments  are  measured  about 
and  the  wing  is  pitching  about  the  leading  edge. 

In  order  to  illustrate  the  result  given  by  equation  (37),  the 
damping  moment  due  to  steady  pitching  has  been  computed  for  a family  of 
airfoils  having  symmetrical  biconvex  parabolic  arc  sections.  The 
results  are  presented  in  figure  6.  The  dashed  line  in  figure  6 repre- 
sents the  theoretical  locus  of  Mach  number  and  leading-edge  angle  for 
bow -wave  detachment  (ref.  24).  Since  the  requirements  of  the  second- 
order  theory  are  violated  to  the  left  of  this  line,  the  damping  curves 
there  are  not  rigorously  correct,  and  therefore  are  terminated  a short 
distance  beyond  the  line. 

To  determine  the  effect  of  thickness  on  the  parameter,  Cm£' , it  is 
assumed  that  the  shape  of  the  indicial  pitching-moment  variation  is  not 
significantly  affected  by  thickness.  When  this  assumption  is  made,  it 
becomes  necessary  only  to  correct  the  initial  and  final  ordinates  of  the 
indicial  curve  Cn^' (o)  and  Cn^' (oa)  and  the  number  of  half-chord  lengths 
traveled  to  reach  steady  state,  aa.  The  response  curve  is  adjusted  to 
fit  these  new  end  values,  so  that  the  area  representing  the  cor- 
rected Cjjj^'  can  then  be  determined. 

The  effect  of  thickness  on  the  final  ordinate  of  the  indicial  curve, 
Cn^' (oa),  has  already  been  determined  (eq.  ( 34) ) . Now  consider  the 
starting  lift  and  moment. 

Assume  first  a thin  flat  wing  to  be  flying  at  zero  angle  of  attack 
in  a uniform  stream  of  density  p0.  At  time  zero,  it  starts  to  sink 
with  downward  velocity  V0a.  The  impulsive  start  causes  a plane  com- 
pression wave  to  be  emitted  from  the  lower  surface  at  the  speed  of  sound, 
which  is  constant  throughout  the  stream.  (The  pressure  and  therefore 
the  temperature  is  sensibly  constant.)  At  the  same  instant,  a plane 
expansion  wave  is  emitted  from  the  upper  surface  and  travels  at  the  same 
velocity  in  the  opposite  direction.  At  the  end  of  an  infinitesimal  time 
period  At,  the  wing  has  sunk  a distance  VoaAt.  The  sound  waves  have 
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moved  a distance  aoAt  in  both  directions.  The  mass  of  fluid  affected 
is  2p0c0aoAt.  Then,  from  Newton's  second  law,  one  may  write_ 

(fu-fz)  At  = 2p0c0aoAt  V0a 

where  fu  and  f i are  the  forces  per  unit  span  on  the  upper  and  lower 
surfaces,  respectively.  Converting  to  coefficient  form,  it  then  appears 
that 


CL(o)  _ fu-fZ  _ ± 

a \ PoVo2coa 

For  the  wing  of  finite  thickness  for  which  the  pressure  varies  along  the 
chord,  a different  result  is  expected.  The  wing  is  assumed  to  have  been 
flying  at  zero  angle  of  attack  for  a period  of  time  long  enough  for  the 
flow  to  have  established  itself  around  the  surface.  In  this  case,  the 
density,  temperature,  and  therefore  the  speed  of  sound,  vary  along  the 
chord  due  to  the  curvature  of  the  surface.  At  time  zero  the  wing  again 
starts  to  sink,  but  the  compression  and  expansion  waves  are  propagated 
at  different  rates  at  various  chordwise  positions  so  that 


At  [fu-fl](x)  = 2 (x)  ^ (x)  PcPoAt  V0a 


and  the  starting  lift  becomes 

PL.(.°i  sJL_JC°-£.  a.dx 
a ^0CO  0 Po  ao 


(38) 


From  reference  24,  the  variation  of  p/pQ  a/aQ  with  Cp  is  found  to  be 


_P_  _a_ 

Po  ao 


(39) 


where  Cp  is  the  pressure  coefficient  at  a point  on  the  surface  when 
the  airfoil  is  at  zero  angle  of  attack.  Again,  from  reference  24, 


C 


P 


P-Po 

«o 


= C 


xou  + C2c 


u 


where  p and  p0  are,  respectively,  the  pressure  at  a point  on  the  sur- 
face of  the  airfoil  and  the  free-stream  static  pressure. 
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For  the  symmetrical  biconvex  parabolic  arc  airfoil  considered 
earlier,  the  local  slope  ou  is  given  by 


o 


u 


25 


2x 
Co  . 


w 


where  5 is  the  thickness  ratio. 

Then  expanding  equation  (39)  in  & binomial  series,  retaining  terms 
through  the  second  order  in  ou,  and  performing  the  integration  indicated 
by  equation  ( 38) , there  results  for  the  starting  lift  coefficient 

cLa,(°)  = ~ |l+0.0^Mo262 [20C2-Mo2C12 ] j (42) 


A similar  procedure  can  be  carried  out  for  the  starting  pitching-moment 
coefficient,  measured  about  the  leading  edge 


Cma 


'(o) 


fC°  P a 
MqCo2  0 Po  “ 


Cma*  (o)  = - ^ |l+0.04Mo2[20C262-Mo2C1252-10C16]  j (43) 

It  is  necessary  also  to  consider  the  effect  of  thickness  on  the  time 
required  for  the  loading  on  the  airfoil  subjected  to  a step  change  in 
angle  of  attack  to  reach  steady  state.  As  was  seen  in  sketches  (n) 
and  (p),  the  time  to  reach  steady  state  is  exactly  equivalent  to  the 
time  required  for  the  wing  to  escape  the  influence  of  the  sound  waves 
emitted  from  the  leading  edge  at  the  start  of  the  motion.  For  purposes 
of  analysis,  the  situation  is  reversed  by  requiring  that  the  wing  remain 
stationary  in  a fluid  moving  with  supersonic  velocity  and  causing 

a sound  wave  to  be  emitted  at  the  leading  edge  of  the  wing.  To  an 
observer  standing  on  the  upper  surface  of  the  airfoil,  there  will  appear 
to  be  two  wave  fronts  whose  velocities  tangent  to  the  surface  are 
[V+a]  and  [V-a],  where  V and  a are  the  local  stream  velocity  and  local 
speed  of  sound,  respectively.  Then  the  time  required  for  the  slower 
wave  front  to  clear  the  airfoil  is  given  by 


^a 


dx 

[V-a]cos[tan“xou] 


Mq  r00  /i+ou2  a* 

vo  Jo  [V/a0-a/a0] 


(44) 
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It  is  therefore  necessary  to  obtain  v/a0  and  a/a0  as  functions  of 
chordwise  location  x.  Again,  in  reference  24  these  quantities  are 
found  to  be  functions  of  the  local  pressure  coefficient  Cp. 


a 


1 + 


0*5) 


Expanding  equations  (45)  through  the  second  power  of 
through  the  second  power  of  au,  there  is  obtained 


and 


n/'l+c 


2 

u 


ta  = 


Mo 

VMo 


1 + - Oi 


_ r l 2 j 

-1  ] JQ  l-K1Cp+K2C 


dx 


(U6) 


where 


O.^OMp  + O.lOMp2 
Mq  - 1 


K 


2 


0.125Mo[Mog-l3  + 0.03Mq4 
Mo  - 1 


The  pressure  coefficient  Cp  is  now  assumed  to  be  that  existing  at 
a point  on  the  surface  when  the  airfoil  is  at  zero  angle  of  attack.  This 
assumption,  of  course,  introduces  some  error  since  the  angle  of  attack 
has  changed  discontinuously  from  zero  to  a at  t = 0.  However,  if  we 
require  that  a be  less  than  &2,  then  the  error  introduced  through 
neglect  of  the  effect  of  the  change  in  a on  the  pressure  coefficient 
is  negligible.  Then,  as  in  the  starting  lift  problem. 


Cp  = Cl0u  + C2aua  (47) 

Inserting  equation  (47)  into  equation  (46),  dividing  numerator  by  denom- 
inator, and  again  retaining  terms  through  the  second  power  in  ou,  there 
results 


Mq 

’a  ' v0  [Mq-i  J 


l+C^k^Oy 


k1C2+[KI2-K2]C12  + 


dx  (48) 
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For  the  symmetrical  biconvex  parabolic  arc  airfoil,  cu  is  replaced  by 
equation  (41),  whereupon  for  this  airfoil  the  time  to  reach  steady- 
state  becomes 


coMo 

V0[Mo-l) 


(1+1  ] 


(49) 


where 


i j*1Cs+[«12-*2]Cl2+ij 

Alternatively,  in  terms  of  the  number  of  half-chord  lengths  of  travel, 

2V0  2Mq  r . , 

oa  = ^ ta  = -J2-  [1+|]  (50) 

c0  Mq-1 

The  effect  of  thickness  on  the  damping  parameter  Cm^'  can  now  be 
estimated  approximately  by  adjusting  the  indicial  pitching-moment  curve 
obtained  from  the  linearized  theory  to  the  end  values  given  by  equa- 
tions (34),  (43),  and  (50),  and  then  measuring  the  area  corresponding  to 
the  "corrected”  Cm^.  The  process  is  indicated  in  sketch  (dd)  below, 
and  the  adjustment  formulas  are  given  by  equations  (51)* 
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<Pl 


Cmao,(<P0^  " ^“a,'  (°)  " ^ 


(og,)  " ACm^'  (o) 


JP fl 

°ar 


(51) 


OaA 


where 


ACma'lo)  = C^' (o)  - C^'Co) 


ACn^' (oa)  = Cm^'  (°a0)  “ Cmc^'^ai) 


and  the  subscripts  o and  l indicate  values  corresponding  to  the  linear- 
ized theory  and  the  adjusted  theory,  respectively.  Notice  also  in 
sketch  (dd)  that  the  area  corresponding  to  the  adjusted  value  of  Cm^' 
is  given  by  the  following  expression: 

Cm^'  = Cm^'  + " °a^  ACn^* (eg)  (52) 


The  parameter  Cm^’  for  pitching  about  the  leading  edge  has  been 
computed,  using  equations  (35)  to  (52)  for  the  same  family  of  airfoils 
as  that  discusse’d  previously  and  the  results  added  to  those  obtained 
for  Cjn^’  (eq.  (37))*  The  total  damping  moment  Cmq*  + cm<x'  is  shown 
plotted  as  a function  of  Mach  number  in  figure  7*  It  is  evident  from 
examination  of  figure  7 that  the  effects  of  thickness  are  small  and  are 
in  the  destabilizing  direction  with  increasing  thickness  ratio.  It 
should  be  noted  that  this  result  does  not  agree  with  those  of  Jones  and 
of  Wylly  in  references  25  and  26  which  indicate  that  the  effect  of  thick- 
ness is  extremely  stabilizing.  In  reference  27,  however,  Van  Dyke  takes 
issue  with  the  results  of  references  25  and  26,  and  offers  an  alternate 
solution  correct  to  the  second  order  in  thickness.  From  this  result, 
the  following  relationship  for  the  damping-moment  coefficient  for  pitch- 
ing about  the  leading  edge  can  be  deduced: 


Cmq'+Cma  - Cn^'+Cm^'  1 + { 


f X3[Mo2N-2] 


2-Mq‘ 


>}] 


(53) 


where 


Cmq'+Cm<x  is  the  damping -moment  coefficient  obtained  from  the 

- Jrt 


linearized  theory  and 
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X = 1 for  biconvex  airfoils 

3 

X = for  double-wedge  airfoils 

The  trend  of  the  damping-moment  coefficient  with  thickness  obtained  from 
this  result  agrees  closely  with  that  obtained  in  the  present  paper. 

EXPERIMENT 
Wind  Tunnel 


The  experimental  investigation  of  the  damping-in-pitch  character- 
istics of  the  various  model  configurations  was  conducted  in  the 
Ames  6-  by  6-foot  supersonic  wind  tunnel.  In  this  tunnel,  the  Mach  num- 
ber can  be  varied  continuously  within  the  subsonic  speed  limits  0.60  to 
choking,  and  at  supersonic  speeds  from  1.20  to  2.00.  The  total  pressure 
can  also  be  changed  at  any  time  within  the  limits  3 pounds  per  square 
inch  absolute  to  atmospheric.  A complete  description  of  the  wind  tunnel 
is  given  in  reference  28. 


Models 


The  pertinent  dimensions  of  the  five  wing-body  configurations 
investigated  in  the  present  report  are  shown  in  figure  8.  All  the  wings 
had  sections  which  were  symmetrical  in  streamwise  planes  and  3 percent 
thick.  Wings  having  leading-edge  sweep  angles  less  than  or  equal  to  45° 
(wings  c,  d,  and  e in  fig.  8)  had  biconvex  circular-arc  sections  in 
streamwise  planes  with  maximum  thickness  at  50-percent  chord.  Wings 
having  leading-edge  sweep  angles  greater  than  45°  (wings  a and  b in 
fig.  8)  had  NACA  0003-63  sections  in  streamwise  planes.  The  wings  were 
identical  in  plan  form  and  section  shape  to  those  investigated  in  the 
series  of  force  tests  reported  in  references  29  to  33*  Strength  consid- 
erations required  that  wings  c,  d,  and  e be  built  of  steel.  Wings  a 
and  b were  built  of  aluminum. 

The  models  were  fitted  with  bodies  such  that  the  distance  from  the 
apex  of  the  body  to  the  wing-body  intersection  was  the  same  for  each 
model.  Also,  the  ratio  of  wing  area  to  body  maximum  cross-sectional  area 
was  17*9  for  each  model.  Forward  of  the  point  of  maximum  radius,  the 
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bodies  were  of  thin  laminated  wood  construction  and  were  identical  in 
shape  to  those  of  references  29  to  33*  The  aluminum  afterbodies  were 
cylindrical  in  shape  and  terminated  at  the  trailing  edges  of  the  wings 
to  permit  the  models  to  be  deflected  with  sufficient  amplitude'  when 
mounted  on  the  support  structure. 


Model  Support  System 


The  support  system  used  in  the  present  investigation  was  the  same 
as  that  developed  for  the  single-degree-of-freedom  free-oscillation 
experiments  of  reference  10.  The  interested  reader  is  referred  to  that 
report  for  a detailed  description  of  the  oscillation  mechanism.  For 
the  present  tests,  two  changes  were  made  in  the  system  which  served  to 
improve  the  accuracy  of  the  data.  The  changes  and  the  reasons  for  them 
are  as  follows: 

1.  The  vertical  flexure  pivots  which  provided  vertical  restraint 
in  the  tests  of  reference  10  were  replaced  by  bearings  after  it  was 
found  that  the  heavier  models  of  this  investigation  caused  the  pivots 
to  twist  laterally  when  the  model  was  given  its  initial  displacement. 

An  undesirable  yawing  oscillation  was  thereby  induced  which  destroyed 
the  linearity  of  the  system.  Installation. of  the  bearings  entirely  elim- 
inated this  yawing  tendency.  The  damping  due  to  friction  in  the  bearings 
was  somewhat  larger  than  that  of  the  flexure  pivots,  but  by  frequent 
checking  of  the  wind-off  damping  and  regular  replacement  of  the  bearings 
it  was  possible  to  maintain  the  tare  damping  well  within  acceptable 
limits . 


2.  A stiffening  strut  was  installed  between  the  sting  support  and 
the  tunnel  ceiling  in  order  to  remove  the  possibility  of  coupling  between 
the  model  oscillation  and  the  resonant  mode  of  the  sting  support.  The 
influence  of  the  strut  on  the  aerodynamic  forces  at  subsonic  speeds  was 
investigated  by  removing  the  strut  and  recording  data  for  one  of  the 
models  at  a frequency  sufficiently  below  the  sting  resonant  frequency  to 
avoid  excessive  sting  vibrations.  Results  obtained  with  the  strut 
installed  and  removed  agreed  within  the  experimental  accuracy. 

By  effectively  eliminating  both  yawing  tendencies  of  the  model  and 
vibrations  of  the  sting  support,  it  can  be  seen  from  comparison  of  the 
present  results  with  those  of  reference  10  that  the  deviations  between  a 
number  of  observations  at  a given  Mach  number  were  markedly  reduced. 
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Scope  of  Tests 


Investigation  of  the  damping-in-pitch  characteristics  of  the  various 
model  configurations  was  conducted  over  a supersonic  Mach  number  range 
of  1.20  to  1.90  and,  where  possible,  over  a subsonic  Mach  number  range 
of  0.60  to  0.90.  The  Reynolds  number  for  the  major  portion  of  the  tests 
was  held  constant  at  1.6  X 10s  per  foot.  Additional  data  for  one  of  the 
models  was  obtained  at  Reynolds  numbers  of  1.0  X 106  per  foot  and 
3.2  x 10®  per  foot. 

All  the  tests  were  conducted  with  the  models  at  a mean  angle  of 
attack  of  0°,  the  angle  being  measured  from  the  mean  line  of  the  sting 
support  to  the  axis  of  the  test  section.  For  each  oscillation,  an  ini- 
tial displacement  of  ±5°  was  imparted  to  the  model  by  means  of  the 
pneumatically  actuated  pawl  arrangement  described  in  reference  10. 

Other  important  variables  are  listed  below  for  the  five  wing-body 
configurations  investigated: 


Model 

Range  of  moment 
of  inertia, 
slug  - ft2 

Axes  of 
rotation, 

% M.A.C. 

Range  of  reduced 
frequency 

A = 2 triangular 

0.0534  - O.O589 

35,  45 

0.012  - 

0.037 

A = 3 triangular 

.0368  - .0402 

25,  35 

.018  - 

.050 

A = 4 triangular 

.0491  - .0512 

25,  35,  45 

.007  - 

.o4i 

A = 3 swept 

.0549  - .0576 

25,  35 

.011  - 

.025 

A = 3 unswept 

.0422  - .0488 

20,  35,  40 

.007  - 

.025 

Reduction  of 

Data 

The  technique  used  in  this  investigation  for  reducing  an  experi- 
mental oscillation-decay  record  to  the  damping-in-pitch  coefficient 
Cmn  + Cm^  was  identical  to  that  described  in  reference  10.  Briefly, 
this  technique  involves  plotting  the  envelope  of  the  oscillation-decay 
curve  against  time  on  semilogarithmic  graph  paper.  If  the  damping  of 
the  system  is  a linear  function  of  the  angular  velocity,  this  plot  will 
be  a straight  line,  the  slope  of  which  is  proportional  to  the  damping 
term.  Calibrations  for  the  model  moment  of  inertia  and  the  damping  due 
to  internal  friction,  and  standard  measurements  for  the  density  and 
velocity  of  the  air  stream  then  enable  one  to  reduce  the  aerodynamic 
damping  term  to  coefficient  form.  The  aerodynamic  restoring  moment  C^ 
can  be  obtained  from  measurements  of  the  frequency  of  oscillation  and 
a calibration  for  the  static  spring  constant  of  the  system.  Results  for 
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the  restoring  moment  about  two  axes  of  rotation  a known  distance  apart 
then  permit  calculation  of  the  lift-curve  slope  CLa» 


Corrections  to  Data 


No  wind-tunnel- wall  corrections  were  made  to  the  subsonic  results 
for  any  of  the  measured  stability  derivatives.  To  the  author's  knowl- 
edge, no  corrections  have  been  developed  for  application  to  the  forces 
acting  on  a finite-span  wing  oscillating  in  a compressible  fluid. 
Tunnel-resonant-frequency  effects  were  investigated  by  use  of  refer- 
ence and  it  was  found  that  the  range  of  frequencies  used  in  this 
investigation  was  well  below  the  lowest  calculated  tunnel  resonant  fre- 
quency. In  view  of  the  fact  that  the  subsonic  damping  moments  measured 
for  the  A = k triangular  wing  in  the  present  investigation  agreed 
reasonably  well  with  the  results  for  a similar  wing  obtained  during  an 
investigation  in  the  Ames  12-foot  wind  tunnel  (ref.  35) t it  is  believed 
that  the  effect  of  the  tunnel  walls  on  the  damping  in  pitch  was  not 
significant. 

A subsonic  correction  could  have  been  made  to  the  static  parameters 
C]^  and  Cjd^.  However,  since  this  correction  would  have  been  very  small, 

and  in  view  of  the  uncertainty  involved  in  applying  static  corrections 
to  the  results  of  dynamic  measurements,  it  was  decided  not  to  make  the 
correction. 

The  effects  of  constriction  of  the  flow  at  subsonic  speeds  by  the 
tunnel  walls  were  taken  into  account  by  the  method  of  reference  36. 

This  correction  amounted  to,  at  most,  a 2-percent  increase  in  the  Mach 
number  and  in  the  dynamic  pressure  over  that  determined  from  a calibra- 
tion of  the  wind  tunnel  without  a model  in  place. 

For  the  tests  at  supersonic  speeds,  no  corrections  were  required  to 
be  made  on  either  the  aerodynamic  measurements  or  the  air-stream 
properties. 

The  reader  is  referred  to  reference  10  for  a discussion  of  the 
precision  of  the  data.  Since  the  method  of  reduction  of  the  data  was 
the  same  as  that  of  reference  10,  the  uncertainty  in  the  recorded  value 
of  a given  damping  coefficient  remained  of  the  order  of  ±0.02.  However, 
by  virtue  of  the  elimination  of  yawing  tendencies  in  the  model  and  sting- 
support  vibrations,  the  standard  deviation  of  a number  of  observations 
at  a given  Mach  number  was  reduced  from  ±0.06  to  ±0*03* 
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RESULTS 

Damping  Coefficients 


The  basic  experimental  results  for  the  damping  coefficient  Cm^+Cm^ 
for  the  five  wing-body  configurations  investigated  are  shown  in  fig- 
ures 9 to  13*  In  cases  where  data  could  be  obtained  at  subsonic  as  well 
as  supersonic  speeds,  these  results  are  shown  plotted  through  the  Mach 
number  range  0.60  to  1.90  for  ready  comparison  of  the  magnitudes  of  the 
coefficients  in  the  two  speed  ranges.  Also  shown  plotted  in  figures  9 
to  13  are  the  theoretical  results  for  the  damping  coefficient  Cm^  + Cn^ 
at  supersonic  speeds  for  the  wings  alone.  In  addition,  theoretical 
results  for  the  triangular  wing-body  combinations,  considering  inter- 
ference effects,  are  shown  in  figures  9 to  11.  The  sources  of  these 
theoretical  results  are  references  11,  16  through  23,  and  37*  The  exper- 
imental results  for  subsonic  speeds  are  not  compared  with  theory,  since 
a rigorous  linearized  theory  for  triangular,  swept,  or  unswept  wings 
oscillating  in  a subsonic  compressible  flow  has  not  yet  been  developed. 

A very  approximate  analysis  for  Cmq  and  Cm^  for  the  aspect  ratio  4 
triangular  wing  was  presented  in  reference  10.  In  this  development,  the 
parameter  was  approximated  using  the  incompressible  theory  of  ref- 

erence 8 for  the  elliptic  wing.  However,  a more  rigorous  analysis  for 
the  elliptic  wing  in  compressible  flow  (ref.  11)  has  indicated  that  the 
result  of  reference  10  for  Cn£  may  be  incorrect.  The  good  agreement 
between  the  experimental  subsonic  results  of  reference  10  and  the  sub- 
sonic theory  of  reference  10  may  be  fortuitous. 


The  results  for  the  three  triangular  wings  (figs.  9 to  ll)  and  the 
swept  wing  (fig.  12)  show  that  in  the  supersonic  speed  range,  in  every 
case,  the  linearized  first-order-in -frequency  theory  provides  a reliable 
guide  for  obtaining  both  the  magnitude  and  trend  with  Mach  number  of  the 
damping  coefficients.  Of  particular  importance  are  the  theoretical 
predictions  of  ranges  of  Mach  number  over  which  the  wings  could  experi- 
ence negatively  damped  oscillations.  The  experimental  results  for 
the  A * 3 triangular  wing  (fig.  10)  are  inconclusive  in  this  regard 
since  the  theoretically  predicted  unstable  range  occurs  in  the  low  super- 
sonic Mach  number  range,  1.0  to  about  1.08,  where  it  was  not  possible  to 
obtain  data.  The  predictions  by  the  theory  of  larger  unstable  ranges 
for  the  A » 4 triangular  wing  and  the  A = 3 swept  wing  were  borne  out 
by  the  experimental  results,  as  shown  in  figures  11  and  12  by  the  exper- 
imented points  plotted  below  the  abcissas.  In  the  absence  of  available 
theoretical  results  for  unswept  wings,  the  experimented  results  for  the 
unswept  wing  are  compared  with  theoretlced  damping  coefficients  for  a 
rectangular  wing  of  aspect  ratio  3 (fig*  13)*  It  is  seen  from  examina- 
tion of  figure  13  that  the  theoretlced  prediction  of  a rather  large 
region  of  instability  for  an  axis  of  rotation  at  0.205  (fig.  13(a))  is 
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borne  out  by  the  experimental  results  but  that  the  theory  fails  to  pre- 
dict the  reversal  in  trend  of  the  damping  coefficients  at  low  supersonic 
Mach  numbers  for  axes  of  rotation  closer  to  the  aerodynamic  center 
(figs.  13(b)  and  (c)).  It  is  clear  from  these  results  that  the  effects 
of  taper  can  be  important  and  should  be  considered  for  an  accurate  theo- 
retical appraisal  of  the  damping- in-pitch  characteristics  of  unswept 
wings. 

It  should  be  mentioned  that  at  a Mach  number  of  0.90,  the  damping 
of  the  A = 4 triangular  wing  (fig.  11)  and  the  A = 3 triangular  and 
swept  wings  (figs.  10  and  12)  were  highly  nonlinear  with  angle  of  attack, 
being  stabilizing  for  angles  of  attack  of  5°  and  decreasing  steadily  to 
zero  at  about  1°.  The  dashed  lines  for  this  Mach  number  shown  in  fig- 
ures 10,  11,  and  12  are  intended  to  indicate  the  range  of  this  nonlinear 
variation.  Further  discussion  of  the  phenomenon  is  withheld  to  a later 
section  of  this  report.  The  above-mentioned  difficulty  was  not  encoun- 
tered with  the  A = 2 triangular  wing,  and  the  value  of  Cm(j  + 
at  Mq  = 0.9  for  this  wing  (fig.  9)  represents  the  damping  throughout 
the  amplitude  range.  Subsonic  results  for  the  unswept  wing  could  not  be 
obtained  due  to  the  highly  erratic  behavior  of  the  aerodynamic  restoring 
moment  which  generally  was  so  strongly  destabilizing  as  to  counterbalance 
the  spring  restoring  moment  and  force  the  model  against  its  stops. 


Reynolds  Number  Effects 


For  the  A = 2 triangular  wing  at  supersonic  speeds  (fig.  9), 
damping  coefficients  have  been  obtained  at  three  Reynolds  numbers  rang- 
ing from  1.18  X 10s  to  3.77  X 10®,  while  at  subsonic  speeds  (fig.  9)  they 
have  been  obtained  at  1.18  X 10®  and  I.89  X 10®.  In  both  cases,  it  is 
evident  that  at  least  within  the  range  of  Reynolds  numbers  at  which  data 
were  recorded  there  is  no  significant  effect  of  scale.  It  is  reasonable 
to  assume  that  a similar  lack  of  scale  effects  exists  for  the  other  wings 
as  well. 


Aeroelaetic  Effects 


In  figures  14  to  23,  the  static  parameters  Cn^  and  CLa  for  the 
five  wings  are  compared  with  force-test  measurements  of  similar  wings 
(refs.  29  to  33)*  Such  a comparison  is  useful  as  an  indication  of  the 
importance  of  aeroelastic  effects  on  the  aerodynamic  properties  of  the 
less  rigid  models  of  the  present  tests.  The  fact  that  at  supersonic 
speeds  there  is  excellent  agreement  between  the  two  experiments  for  the 
lift-curve  slope  C^»  whereas  the  present  results  for  are 
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consistently  smaller, is  attributed  primarily  to  the  difference  in  the 
body  shapes  of  the  models  of  the  two  investigations.  In  the  force-test 
investigation,  the  afterbodies  were  carried  beyond  the  trailing  edges  of 
the  wings,  which  would  tend  to  shift  the  centers  of  pressure  of  the  wing- 
body  combinations  farther  rearward  while  leaving  the  total  lift  essen- 
tially unchanged.  From  this  reasoning,  it  is  concluded  that  aeroelastic 
effects  on  the  aerodynamic  coefficients  of  the  present  tests  are  of 
secondary  importance.  At  subsonic  speeds,  there  are  differences  in  mag- 
nitude of  both  Cj^  and  Cn^.  However,  in  view  of  the  fact  that  no  large 
effects  of  aeroelasticity  were  found  at  supersonic  speeds  where  the 
dynamic  pressure  is  greatest,  it  is  believed  that  the  differences 
in  Cj^  and  Cn^  at  subsonic  speeds  are  not  due  to  aeroelastic  effects. 

It  must  therefore  be  stated  that  at  the  present  time  the  differences  at 
subsonic  speeds  are  not  understood. 


Transfer  of  Axes 


Sufficient  information  has  been  given  in  figures  9 to  23  to  permit 
the  transfer  of  the  damping-in-pitch  results  at  supersonic  speeds  to 
axes  of  rotation  other  than  those  used  during  the  tests.  Such  a trans- 
fer, if  valid,  of  course  greatly  increases  the  range  of  applicability 
of  the  given  data.  The  validity  of  the  transfer  equation  (given  in 
appendix  A)  depends  primarily  on  the  invariance  of  the  stability  deriv- 
atives with  changes  in  angle  of  attack  and  angular  velocity  and  could  be 
best  checked  by  experimental  means.  Checks  on  the  method  are  provided 
in  figures  11(c)  and  13(c)  for  the  A = 4 triangular  wing  and  A = 3 
unswept  wing,  respectively.  Shown  in  these  figures  are  the  experimental 
damping  coefficients  obtained  during  the  investigation  and  coefficients 
for  the  same  axis  position  computed,  by  means  of  the  transfer  equation, 
from  the  results  of  experiments  at  two  other  axis  positions.  It  is 
evident  that  the  transferred  results  portray  the  actual  data  with  suffi- 
cient accuracy  to  establish  the  validity  of  the  method. 


DISCUSSION 

Effect  of  Aspect  Ratio 


It  will  be  remembered  from  the  discussion  in  the  theoretical  section 
of  this  report  that  the  linearized  theory  indicates  that  a decrease  in 
aspect  ratio  has  a highly  stabilizing  effect  on  the  damping  in  pitch  of 
a triangular  wing  with  subsonic  leading  edges.  This  indication  is  con- 
firmed in  figure  24,  wherein  the  damping  in  pitch  is  shown  as  a function 
of  Mach  number  for  the  three  triangular  wings  having  their  axes  at  0.35c. 
It  is  noteworthy  that  the  A = 2 wing  shows  no  tendency  toward  dynamic 
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instability  in  the  Mach  number  range  1.2  to  1.9>  that  for  Mq<1.2  the 
trend  of  the  damping  coefficients  of  the  A = 3 wing  is  toward  insta- 
bility, and  that  a range  of  instability  exists  for  the  A = 4 wing.  It 
is  of  further  interest  to  note  that  for  Mach  numbers  greater  than  1.67, 
in  which  range  the  leading  edges  of  the  A = 3 and  A = 4 wings  are 
supersonic,  the  damping  coefficients  for  these  wings  are  essentially  the 
same.  This  serves  as  a partial  confirmation  of  the  theoretical  predic- 
tion of  reference  19,  wherein  it  is  shown  that  the  damping  coefficients 
of  triangular  wings  having  supersonic  leading  edges  are  independent  of 
aspect  ratio. 

A more  comprehensive  way  of  comparing  the  tendencies  of  the  wings 
toward  dynamic  instability  would  be  to  show  the  entire  ranges  of  axis 
positions  and  Mach  numbers  over  which  dynamic  instability  is  possible, 
that  is,  to  plot  the  stability  boundary  curves  for  the  three  wings. 

In  the  Mach  number  range  1.2  to  1.9>  however,  only  the  A = 4-  wing  has 
a region  of  instability  and  this  curve  is  shown,  compared  with  the  theo- 
retical curve,  in  figure  25. 

The  preceding  comparison  (fig.  24-)  was  useful  primarily  for  verify- 
ing the  theoretical  prediction  regarding  the  role  played  by  aspect  ratio 
in  determining  the  tendency  of  a triangular  wing  toward  dynamic  instabil- 
ity. A comparison  of  the  actual  magnitudes  of  the  damping  moments  was 
masked,  however,  by  the  fact  that  the  coefficients  for  the  three  wings 
were  referred  to  their  own  respective  characteristic  lengths.  In  the 
following  comparison,  an  attempt  is  made  to  overcome  this  difficulty  by 
posing  the  question:  Given  three  triangular  wings  of  aspect  ratios  2,  3 > 

and  4,  of  equal  area,  and  required  to  have  the  same  restoring  moment 
(in  ft-lb)  at  a Mach  number  of  1.4,  how  do  the  physical  magnitudes  of 
the  damping  moments  compare?  For  this  purpose,  the  aspect  ratio  3 wing 
is  chosen  as  a standard,  required  to  have  a static  margin  of  0.05 
at  M0  = 1.4,  and  the  damping  moments  of  all  the  wings  referred  to  the 
M.A.C.  of  the  aspect  ratio  3 wing.  The  comparison  on  this  basis  is 
shown  in  figure  26.  It  is  evident  from  examination  of  figure  26  that  on 
this  basis  also,  the  effect  of  decreasing  the  aspect  ratio  is  beneficial 
to  the  damping  in  pitch,  the  aspect  ratio  2 wing  having  the  largest 
damping  moments  throughout  the  range  of  test  Mach  numbers,  and  the  aspect 
ratio  4 wing  the  smallest. 


Effect  of  Plan-Form  Shape 


For  this  comparison,  use  is  made  of  the  results  for  the  aspect 
ratio  3 triangular,  swept,  and  unswept  wings.  First,  the  prediction 
made  in  the  section  entitled  "Theory"  regarding  the  tendencies  of 
three  similar  wings  toward  dynamic  instability  is  examined  by  comparing 
their  stability  boundaries.  From  the  preceding  section,  it  will  be 
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recalled,  however,  that  no  region  of  instability  exists  for  the  A = 3 
triangular  wing  in  the  Mach  number  range  1.2  to  1.9*  It  is  apparent, 
therefore,  that  as  predicted  by  the  theory  the  triangular  wing  is 
superior  in  this  regard.  The  stability  boundary  curves  for  the  swept 
and  unswept  wings,  with  axis  position  referred  to  the  leading  edge  of 
the  M.A.C.,  are  compared  in  figure  27.  Again,  the  theoretical  predic- 
tion is  borne  out  for,  as  seen  in  figure  27,  a considerably  larger  region 
of  instability  exists  for  the  unswept  wing. 

Next,  the  magnitudes  of  the  damping  moments  are  compared  by  choosing 
the  A ■ 3 triangular  wing  as  a standard,  with  a static  margin  of  0.05 
at  Mo  = 1.4,  and,  as  in  the  previous  section,  requiring  that  the  wings 
have  equal  restoring  moments  at  this  Mach  number.  The  comparison  is 
shown  in  figure  28.  On  this  basis,  it  is  again  apparent  that  the  trian- 
gular wing  is  the  superior  wing,  followed  in  order  by  the  swept  and 
unswept  wings. 

Notice  in  figures  2 6 and  28  that  since  the  damping  coefficients  of 
all  five  wings  investigated  are  referred  to  the  M.A.C.  of  the  A = 3 
triangular  wing  and  since  all  wings  have  equal  restoring  moments  at 
Mo  = 1.4,  the  results  of  figures  26  and  28  may  be  compared  directly. 

Here  it  appears  that  the  A » 2 triangular  wing  is  the  most  desirable 
wing  from  a longitudinal -dynamic -stability  standpoint,  and  the  A = 3 
unswept  wing  the  least  desirable.  As  mentioned  previously,  however, 
it  should  not  be  inferred  that  the  result  of  this  comparison  implies 
a recommendation  for  the  use  of  triangular  rather  than  unswept  wings. 
Obviously,  the  addition  of  tail  surfaces  could  alter  considerably  the 
relative  damping-in-pitch  merits  of  the  wings  investigated. 


Effect  of  Thickness 


It  has  been  shown  by  the  results  of  experiments  (figs.  9 to  12) 
that  the  supersonic  linearized  potential  theory  provides  a reliable 
basis  for  predicting  the  low-frequency  damping-in-pitch  characteristics 
of  thin  wings  of  finite  span.  Considerable  doubt  has  been  shed,  however, 
over  the  ability  of  the  linearized  theory  to  predict  the  dynamic  behavior 
of  wings  that  cannot  be  classified  as  thin.  In  recent  years,  a number 
of  reports  have  been  issued  (refs.  25,  26,  and  38)  which  indicate  that 
second-order  thickness  effects  have  a profoundly  stabilizing  effect  on 
the  damping  in  pitch  of  two-dimensional  wings.  In  particular,  refer- 
ence 26  indicates  that  by  increasing  the  thickness  of  an  infinite-span 
wing  from  zero  to  only  4-1/2  percent,  the  rather  large  region  of  insta- 
bility predicted  by  the  linearized  theory  is  completely  eliminated.  In 
contrast  to  these  results,  the  approximate  analysis  of  this  paper 
(see  Theory)  Indicates  that  thickness  effects  are  relatively  small,  the 
destabilizing  effect  on  Cn^  being  partially  canceled  by  a stabilizing 
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effect  on  Cm^.  In  view  of  these  conflicting  results,  it  is  of  some 
interest  to  obtain  additional  experimental  data  which  would  serve  to 
clarify  the  issue.  Of  particular  interest  are  the  effects  of  thickness 
for  finite-span  wings.  In  this  case,  a limited  comparison  can  be  made 
between  the  experimental  results  for  the  3-percent  A = 4 triangular 
wing  of  this  report  and  the  results  of  reference  10  for  two  6-percent 
triangular  wings  having  the  same  aspect  ratio.  Aside  from  the  differ- 
ences in  profile  and  small  differences  in  the  dimensions  of  the  bodies, 
the  models  and  test  conditions  of  the  two  experiments  were  almost  iden- 
tical. The  comparison  is  shown  in  figure  29.  It  is  seen  that  for  this 
wing,  the  effects  of  thickness  are  small  and  are  in  the  direction  of 
decreasing  stability  with  increasing  thickness.  Of  primary  importance 
is  the  fact  that  all  the  wings  exhibit  regions  of  Mach  number  over  which 
dynamic  instability  is  experienced. 

Although  admittedly  the  experimental  evidence  in  figure  29  is  not 
conclusive,  nevertheless  in  the  light  of  these  results  it  is  difficult 
to  accept  those  results  which  indicate  that  large  stabilizing  effects  of 
thickness  are  present  on  infinite-span  wings,  especially  in  view  of  the 
fact  that  the  behavior  of  a triangular  wing  with  supersonic  edges  is  in 
many  ways  quite  similar  to  that  of  an  infinite-span  wing. 


Effect  of  Nonlinearities 


From  examination  of  the  experimental  oscillation-decay  records  at 
both  low  subsonic  and  all  supersonic  speeds,  it  was  evident  that  the 
assumption  of  a linear  second-order  system  was  well  justified.  Aside 
from  a moderate  amount  of  scatter,  a decay  curve  could  be  fitted  by  an 
exponential  curve  over  a wide  range  of  amplitudes,  verifying  that  the 
damping  parameter  was  essentially  independent  of  amplitude.  At  high  sub- 
sonic speeds,  however,  this  linear  behavior  was  no  longer  true  for  most 
of  the  wings  tested.  As  an  example,  figure  30  shows  the  oscillation- 
decay  record  for  the  A = 3 triangular  wing  at  Mq  = 0.90.  It  is  seen 
that  after  an  initial  displacement  the  amplitude  of  oscillation  quickly 
dies  out  to  a low  level,  but  there  sustains  itself  indefinitely. 

Similar  phenomena  have  also  been  reported  in  the  forced  oscillation 
experiments  of  reference  35  for  a wing-body  combination  having  a trian- 
gular wing  of  aspect  ratio  4 and  in  flight  tests  of  various  high-speed 
aircraft,  particularly  tailless  aircraft.  The  deleterious  effects  of 
such  sustained  oscillations  on  the  qualities  of  an  aircraft  as  a gun 
platform  or  from  a structural  standpoint  are  obvious.  It  is  therefore 
of  considerable  importance  to  gain  an  understanding  of  the  phenomenon, 
with  a view  toward  advancing  means  of  eliminating  it. 

At  the  present  time,  insufficient  information,  either  theoretical 
or  experimental,  is  available  to  enable  a complete  and  authoritative 
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description  of  the  physical  event;  a tentative  hypothesis  is  therefore 
advanced  instead,  which  would  appear  to  fit  the  observed  facts.  It  must 
be  emphasized,  however,  that  in  the  light  of  additional  information,  the 
assumptions  and  conclusions  drawn  here  may  require  considerable  revision. 

Specifically,  the  observed  facts  are  these:  (l)  The  damping 

coefficient  is  highly  nonlinear,  being  negative  or  stabilizing  for  angles 
of  attack  greater  than  about  1°,  and  zero  at  1°.  The  possibility  that 
the  nonlinearity  is  caused  entirely  by  scale  effects  and/or  wind-tunnel- 
wall  interference  is  ruled  out  by  the  fact  that  similar  behavior  is 
observed  in  full-scale  flight  tests.  (2)  The  frequency  of  oscillation 
is  essentially  constant  throughout  the  amplitude  range,  implying  that 
the  restoring  moment  is  a linear  or  nearly  linear  function  of  the  angle 
of  attack. 

From  these  observations,  it  is  hypothesized  that  the  characteristic 
differential  equation  governing  the  motion  is  of  the  form 

la  + <p(a,,&)  a + K2a  = 0 

where  cp(a,a)  represents  the  damping  coefficient  as  a nonlinear  function 
of  the  angle  of  attack  and  angular  velocity,  and  I and  K are  constants. 
It  is  therefore  necessary  to  search  for  a mechanism  that  can  affect  the 
damping  moment  to  a much  larger  degree  than  it  does  the  restoring  moment. 
In  this  regard,  it  is  pertinent  to  first  review  some  of  the  character- 
istics of  wings  in  steady  transonic  flow.  It  is  well  known  that  at  high 
subsonic  Mach  numbers,  the  essential  feature  of  the  flow  is  the  appear- 
ance of  shock  waves  on  the  wing  as  the  speed  of  air  over  the  wing  surface 
exceeds  the  local  speed  of  sound.  These  shock  waves  cause  large  changes 
in  the  pressure  and  also  promote  flow  separation  near  the  shock  wave, 
due  to  the  presence  of  adverse  pressure  gradients.  The  effects  of  these 
disturbances  on  the  pitching-moment  characteristics  of  an  unswept  wing 
of  aspect  ratio  3 are  shown  in  reference  39  and  reproduced  as  curve  (a) 
in  sketch  (ee) . Notice  that 
the  slope  of  the  pitching 
moment  is  positive  for  small 
angles  of  attack,  indicating 
that  the  aircraft  would  be 
statically  unstable  in  that 
range.  It  was  also  illus- 
trated in  reference  39  that 
by  reducing  the  curvature  of 
the  wing  profile  near  the 
trailing  edge,  thereby  Sketch  (ee) 

reducing  the  tendency  toward  separation,  the  positive  slope  through  zero 
could  be  eliminated.  This  result  is  shown  as  curve  (b)  in  sketch  (ee). 
By  virtue  of  the  beneficial  effect  of  reducing  separation,  it  is  argued 
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here  that  in  the  absence  of  separation  and  strong  shocks  on  the  basic 
wing,  the  variation  of  pitching  moment  with  angle  of  attack  would  be 
nearly  linear,  as  shown  by  curve  (a)  in  sketch  (ff),  and  that  the  effects 

of  these  disturbances  can 
be  lumped  into  an  additional 
moment  variation,  curve  (b) , 
such  that  the  sum  of  (a) 
and  (b)  then  gives  the 
observed  nonlinear  variation. 


For  triangular  wings, 
the  reversal  in  slope  of  the 
pitching  moment  at  small 
angles  of  attack  does  not 
appear  (see  refs.  29,  30, 
and  31) • This,  however,  does  not  preclude  the  presence  of  the  nonlinear 
variation  due  to  separation.  For  a triangular  wing  with  center  of  grav- 
ity at,  say,  the  quarter  M.A.C.  point,  the  slope  of  the  pitching-moment 
curve  can  be  quite  large.  As  seen  in  sketch  (gg),  the  linear  varia- 
tion (a)  can  then  greatly  outweigh  the  effect  of  the  nonlinear  curve  (b), 
with  the  result  that  a scarcely  evident  nonlinearity  appears  in  the 
combined  curve. 


Sketch  (gg) 


The  preceding  examples  were  taken  from  the 
results  of  static  measurements,  where  the  angle 
of  attack  is  simply  the  angle  between  the  chord 
line  of  the  wing  and  the  free-stream  direction. 
When  the  wing  is  oscillating,  however,  each 
point  | measured  from  the  axis  of  rotation 
experiences  an  additional  angle  of  attack 
a|/vQ,  due  to  the  angular  velocity.  Then  if 
the  steady-state  pitching-moment  variation 
caused  by  the  separation  (curve  (b)  in 
sketch  (gg))  is  approximated  in  the  range 

0<a<amax  ^y 


Cm(  a)  = aa  - ba3 


in  the  unsteady  case  it  is, 
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where  now  | is  assumed  to  be  the  distance  from  the  axis  to  the  point 
at  which  the  additional  lift  due  to  separation  is  concentrated.®  Now, 
since  for  slow  frequencies  <xt/v0  is  much  smaller  than  a,  all  but 
first-order  terms  in  a are  neglected  so  that  equation  (5M  becomes 

Cm  = a [a-ba2]  + ~ [a-3ba2]  (55) 

vo 

The  characteristic  equation  of  motion  then  becomes 


where  [Cmq+Cm&]  and  Ciuq,  are  the  (constant)  stability  derivatives  which 
would  be  present  in  the  absence  of  separation.  Thus,  it  appears  that 
although  the  effect  of  the  nonlinearity  on  the  restoring  moment  can  be 
insignificant,  the  possibility  still  remains  that  the  damping  moment  can 
be  profoundly  affected  in  the  event  that  [Cmq+Cma,]  is  sufficiently  small. 
Notice  that  when  [Cmq+Cm&]  + 2£a/c  is  greater  than  zero  and  if  the  non- 
linear term  in  the  restoring  moment  can  be  ignored,  equation  (56)  takes 
on  the  form 


a - 25a(l-|ia2)  + K2a  = 0 (57) 

Equation  (57)  is  then  recognized  as  being  the  well-known  Van  Der  Pol 
equation  of  nonlinear  mechanics.  It  is  evident  that  for  small  values 
of  a the  damping  term  is  negative,  leading  to  a divergent  oscillation, 
whereas  for  larger  a,  the  damping  term  is  positive.  A stable  regime 
therefore  will  exist  where  a = yi 7c  and  oscillations  of  either  large 
or  small  amplitude  will  converge  to  that  regime. 

Thus,  by  the  assumptions  of  an  additional  nonlinear  pitching-moment 
variation  caused  by  flow  separation  and  a sufficiently  small  damping 
moment  in  the  absence  of  separation,  the  observed  phenomenon  of  a small- 
amplitude  self-sustained  oscillation  can  be  hypothetically  explained  as 
being  due  to  the  destabilizing  effect  of  the  nonlinearity  in  the  restor- 
ing moment  on  the  damping  moment. 

X 1 1 ■ ■ ■ ■ ■ ■ I - — 1 ■ ■ 1 . . ■-  ' 

"It  is  recognized  that  usually  the  nonlinear  effects  of  shock-wave 
boundary-layer  interaction  and  flow  separation  are  markedly  reduced 
during  nonsteady  motions.  However,  for  the  low-frequency  oscillations 
of  the  present  tests  (of  the  order  of  one  cycle  per  100  chord  lengths 
of  travel),  it  is  believed  that  the  effects  of  separation,  though  per- 
haps less  severe  than  for  the  steady  case,  are  nevertheless  still 
present. 
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CONCLUSIONS 

The  results  of  a theoretical  and  experimental  investigation  of  the 
single-degree-of-freedom  damping  in  pitch  of  a series  of  low-aspect-ratio 
wing-body  combinations  made  at  subsonic  and  supersonic  speeds  lead  to  the 
following  conclusions: 

1.  In  the  Mach  number  range  1.2  to  1.9*  theoretical  and  experimen- 
tal values  of  the  damping-in-pitch  coefficient  C^  + were  in  good 

agreement  for  triangular  wings  of  aspect  ratios  2,  3*  and  4,  and  swept 
wing  of  aspect  ratio  3* 

2.  Theoretical  predictions  of  the  existence  of  ranges  of  center- 
of-gravity  positions  for  a range  of  Mach  numbers  greater  than  1.2  over 
which  dynamic  instability  may  be  expected  for  the  aspect  ratio  4 trian- 
gular wing  and  the  aspect  ratio  3 swept  and  unswept  wings  were  confirmed 
by  the  experimental  results. 

3.  The  prediction  by  the  theory  of  the  beneficial  effect  on  the 
damping  in  pitch  at  supersonic  speeds  of  a reduction  in  aspect  ratio  was 
borne  out  by  the  results  of  experiments  for  the  triangular  wings  having 
aspect  ratios  2,  3*  and  4.  In  the  Mach  number  range  1.2  to  1.9*  the 
aspect  ratio  2 wing  had  the  largest  damping  moments  throughout  the  range 
of  Mach  numbers  and  the  aspect  ratio  4 wing,  the  smallest. 

4.  Experimental  results  for  the  stability  boundaries  in  the  super- 
sonic speed  range  of  three  wings  of  aspect  ratio  3 having  triangular, 
swept,  and  unswept  plan  forms  confirmed  the  theoretical  prediction 
regarding  the  relative  magnitudes  of  the  region  of  Mach  number  and 
center-of-gravity  position  in  which  dynamic  instability  could  be  experi- 
enced. In  the  Mach  number  range  1.2  to  1.9#  no  region  of  instability 
existed  for  the  triangular  wing.  The  region  of  instability  for  the 
unswept  wing  was  considerably  larger  than  for  the  swept  wing. 

5*  The  effects  of  profile  thickness  on  the  damping  in  pitch  at 
supersonic  speeds  of  triangular  wings  of  aspect  ratio  4 and  thickness 
ratios  of  3 percent  and  6 percent  were  found  to  be  small  and  in  the 
direction  of  decreasing  stability  with  increasing  thickness  ratio. 

6.  The  occurrence  at  a Mach  number  of  0.9  of  small-amplitude  self- 
sustained  oscillations  of  the  triangular  wing  of  aspect  ratio  4 and  the 
triangular  and  swept  wings  of  aspect  ratio  3 was  attributed  to  a desta- 
bilizing effect  on  the  damping  moment  of  nonlinearities  in  the  aero- 
dynamic restoring  moment. 

Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
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APPENDIX  A 
TRANSFER  OF  AXES 


The  problem  to  be  solved  may  be  posed  as:  Given  the  damping 

coefficient  Cm^  + and  the  restoring-moment  coefficient  about 

two  axes  of  rotation,  a known  distance  apart,  find  the  damping  coefficient 
about  an  arbitrary  axis. 


For  single-degree-of -freedom  rotary  oscillations  about  axis  (l)  the 
damping  coefficient  is  written  as 


Cmq  + - fcm^  + Cm^  + 2xx  — 

Ji  L JQ  L JQ  L J0 


2*1*^  (Al) 


where  the  axis  (o)  is  chosen  as  a reference  axis,  £(3^  + Cm^j^is  known, 

and  xx  is  the  known  (nondimensionalized)  distance  of  axis  (l)  ahead  of 
axis  (o).  For  pitching  about  axis  (2),  the  damping  coefficient  is 


where 


[Cmq  + Si]  = [cm<1  + Cm-]  + 2x^0^]  - + C^j  - 2x^0^ 

[cmq  + Si], 


(A2) 


and  X2  are  known. 


*q  ■ ~m& 

The  relationship  of  the  axes  (l)  and  (2)  to  the  reference  axis  (o)  is 
shown  in  sketch  (hh).  Note  that  if  the  reference  axis  (o)  is  chosen  to 
be  coincident  with  the  axis  passing  through 
the  aerodynamic  center,  then  the  quantities 

2xi|  and  2x  J Cn^  in  equations  (Al)  _ 

^ *^o  *-  -*o  //j  (g)  (oj 

and  (A2)  are  zero.  Since  C™  and  Cnu 

““-l  ^2  Sketch  (hh) 

are  known,  the  parameters  Cj^  and  may  be  determined  from 


Si " Si 

Cr  = _ 

^ *1  ~ X2 

S]0'Si+Sl% , 


(A3) 
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determined,  it  is  now  possible  to  solve  equa- 
tions (Al)  and  (A2)  simultaneously. 


With  Cj^  and 


Let 


[°mq  + Cmd 

£°®q  + Cffld 


= a 


= b 


-‘z 


"“a 


= c 


-*o 


Then 


CLq  + CI* 

- — 1 

_ 

o 

i-b  - c[xx  - x2]  + d[xi2  - x22] 


Cmq  + cm^ 


*1  ~ X2 


= a — xx [c-e ] + dxx2  = f 


= e 


Jo 


(A4) 


The  damping  coefficient  about  cm  arbitrary  axis,  with  x referred  to 
axis  (o)  and  measured  as  positive  ahead  of  axis  (o),  is  therefore. 


Cmq  + Cm^ 


_ -2, 

f + xc  - xe  - x d 


(A5) 
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Mach  number,  M0 

Figure  2.-  Comparison  of  theoretical  single -degree -of -freedom 
short -period  pitching  stability  boundaries  at  supersonic 
speeds  for  three  wings  of  aspect  ratio  3,  having 
triangular,  swept,  and  rectangular  plan  forms . 
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Figure  J.-  The  offset  of  froquoncy  on  the  thooroticol  singto - 
dogroo -of-  froodom  short -poriod  pitching  stability 
boundaries  at  supersonic  speeds  for  a triangular  wing 
of  aspect  ratio  4 . 
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Figure  4.  - The  offset  of  frequency  on  tho  thooro float  single- 
degree -of -freedom  short -poriod  pitching  stability 
boundaries  at  supersonic  speeds  for  a rsetengu/ar 
wing  of  aspect  retie  4. 
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Reduced  frequency,  k 

Figure  5.  - The  variation  with  reduced  frequency  of  the 
single  - degree -of-  freedom  rotary  damping  - moment 
coefficient  for  the  two-dimensional  wing  at  Mach 
numbers  O,  0.50,  and  0.60.  Axis  at  leading  edge. 
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Figure  6.-  The  effect  of  thickness  on  the  damping -in -pitch  coefficient  Cm<f  for  a two-dimensional 
wing  having  a symmetrical  biconvex  parabolic- arc  section.  Axis  at  leading  edge. 


Equation  of  fore  bodies : r • r0ll  - (!—  -yv  J 


Figure  8.-  Sketch  of  models  investigated . 
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Figure  9.—  Export marital  damping- in-pitch  coafficiants  for  th»  wing -body 
combination  having  a triangular  wing  of  aspect  ratio  2. 
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(a)  Axis  at  0.25S. 


Figure  10.—  Experimental  damping -in -pitch  coefficients  for  the  wing - 
body  combination  having  a triangular  wing  of  aspect  ratio  3. 
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Mach  number,  Mu 
(a)  Axis  at  0.256. 

Figure  II.-  Experimental  damping- in- pitch  coefficients  for  the  wing - 
body  combination  having  a triangular  wing  of  aspect  ratio  4 
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Figure  if.—  Concluded . 
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Mach  number,  M0 
(a)  Axis  at  0.256. 


Figure  12.  - Experimental  damping- in- pitch  coefficients  for  the  wing - 
body  combination  having  a swept  wing  of  aspect  ratio  3. 
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Mach  number,  M0 
(b)  Axis  at  0.352. 


Figure  !2.  - Concluded. 


CONFIDENTIAL 


88 


CONFIDENTIAL 


NACA  DM  A52L04a 


MM 

nw^^ww*  /wwwP»| 

ft)  Ante  tf  0.2  OS. 


Figure  13.- Experimental  damping -to- pitch  coefficients  for  the  wing  - 
body  combtootion  boring  on  unswept  wing  of  ospoct  ratio  3. 
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(b)  Axis  at  0.352. 


Figure  /J.  — Continued . 
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Figure  14.—  Experimental  pitching  - moment -curve  slopes  for  the  wing 
body  combination  having  a triangular  wing  of  aspect  ratio  2. 
Axes  of  rotation  at  0.35C  and  0.45 E. 
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Mach  number,  Mm 


Figure  15.-  Experimental  pitching  - moment  -curve  slopes  for  the  wing- 
body  combination  having  a triangular  wing  of  aspect  ratio  3. 

Axes  of  rotation  at  0.2  56  and  0.356. 
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Figure  16.- Experimental  pitching -moment  -curve  slopes  for  the  wing- 
body  combination  having  a triangular  wing  of  aspect  ratio  4. 
Axes  of  rotation  at  0.256,  0.356,  and  0.456. 
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Figure  17.  — Experimental  pitching -moment  -curve  slopes  for  the  wing- 
body  combination  having  a swept  wing  of  aspect  ratio  3. 

Axes  of  rotation  at  0. 25 E and  0.35E. 


CONFIDENTIAL 


3X 


NACA  PM  A52L0Ua 


CONFIDENTIAL 


95 


Mach  number,  AC 


Figure  18.—  Experimental  pitching -moment  - curve  s/opes  for  the  wing • 
body  combination  haring  on  unswept  wing  of  aspect  ratio  3. 
Axes  of  rotation  at  0.203,  0.353 , and  0.403, 
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Figure  19.—  Experimental  lift -curve  slopes  for  the  wing -body 
combination  having  a triangular  wing  of  aspect  ratio  2, 
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Figure  20— Experimental  lift  - curve  slopes  for  the  wing -body 
combination  having  a triangular  wing  of  aspect  ratio  3. 
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Figure  2t.  —Experimental  lift  - curve  slopes  for  the  wing-body 
combination  hoeing  a triangular  wing  of  aspect  ratio  4. 
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Figure  22.— Experimental  lift  -curve  slopes  for  the  wing- body 
combination  having  a swept  whig  of  aspect  ratio  3. 
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Figure  23.-  Experimental  lift -curve  slopes  for  the  wing -body 
combination  having  an  un swept  wing  of  aspect  ratio  3. 
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Figure  25. - Comparieon  of  theoretical  and  experimental  single  - 
degree -of- freedom  short-period  pitching  stability 
boundaries  for  a wing- body  combination  hoeing  a 
triangular  wing  of  aspect  ratio  4. 
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Figure  26.  — Comparison  of  magnitudes  of  experimental  rotary  damping- moments  for  three 

wing  - body  combinations  having  triangular  wings  of  aspect  ratios  2,  3,  and  4. 
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Figure  27.-  Comparison  of  experimental  single  - degree  - of - 
freedom  short  - period  pitching  stability  boundaries  for  two 
wing  - body  combinations  having  swept  and  unswept  wings 

of  aspect  ratio  3 . 


CONFIDENTIAL 


CONFIDENTIAL 


Figure  28.  — Comparison  of  magnitudes  of  experimental  rotary  damping  moments  for  three 
wing -body  combinations  haring  triangular , swept,  and  unswept  wings  of  aspect  ratio  3. 
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Figure  29.- Comparison  of  experimental  damping -in •pitch  coefficients  for  three  wing  -body 
combinations  haring  triangular  wings  of  aspect  ratio  4 and  thickness  ratios  of  3 and  6 percent. 


Oscillograph  record  at  Mach  number  0. 90  for  the  w/hg  - body  combination 
having  a triangular  wing  of  aspect  ratio  3.  Axis  at  0.256. 
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